Quantum gravimeters and gradiometers

ABSTRACT

A method of measuring an acceleration using one or more trapped BECs and at least two atomic field modes within those BECs is described. A density distribution of at least one of the one or more BECs is modified by the acceleration, which leads to a change of the field modes&#39; time evolution. The method comprises: selecting two modes that are differently affected by the acceleration; and measuring the acceleration-induced difference; and using the measured acceleration-induced difference to infer the acceleration. Also described is a method of measuring a field gradient using one or more trapped BECs. Each of the trapped BECs has a density distribution, and at least two atomic modes within those BECs, wherein the density distributions of the one or more BECs are modified by the field, which leads to a change of the atomic modes&#39; time evolution. The method comprises selecting two atomic modes that are differently affected by the field; and measuring a field-induced difference between the selected atomic modes; and using the measured field-induced difference to infer a gradient of the field.

FIELD OF THE INVENTION

The present invention relates to the measurement of acceleration using quantum effects, in particular acceleration due to gravity, and particularly but not exclusively to quantum gravimeters comprising Bose-Einstein condensates (BECs). The present invention also or alternatively relates to the measurement of a field gradient using quantum effects, in particular a gravity gradient, and particularly but not exclusively to quantum gradiometers comprising BECs.

BACKGROUND TO THE INVENTION

Quantum gravimeters exploit quantum properties to measure local gravitational fields with high precision. Applications range from finding pipework under London to navigation and search for oil. Industry and science would strongly benefit not only from a precise instrument but also from a portable one.

Current quantum gravimeters use atom interferometers that measure between a few meters and centimetres and cannot be further miniaturised without losing precision. This is because their precision depends on the interferometer arm length through the atomic time of flight in the interferometer.

Gravitational waves are postulated in the theory of general relativity as distortions in the local gravitational field. It has been shown that such gravitational waves would produce collective excitations (called phonons) in a BEC, and that those phonons are detectable under appropriate conditions.

In C. Sabin, D. E. Bruschi, M. Ahmadi, and I. Fuentes, New J. Phys 16 085003 (2014) it is shown that gravitational waves create phonons in a Bose-Einstein condensate (BEC). A traveling spacetime distortion produces particle creation resonances that correspond to the dynamical Casimir effect in a BEC phononic field contained in a cavity-type trap. It is proposed to use this effect to detect gravitational waves. The amplitude of the wave can be estimated by applying recently developed relativistic quantum metrology techniques. The optimal precision bound on the estimation of the wave's amplitude is provided. Finally, it is shown that the parameter regime required to detect gravitational waves with this technique could be, in principle, within experimental reach in a medium-term timescale.

M. Ahmadi, D. E. Bruschi, C. Sabin, G. Adesso and I. Fuentes, Sci. Rep. 4, 4996 (2014) presents a framework for relativistic quantum metrology, in which techniques for the application of metrology to quantum field theory are introduced. This framework allows for high precision estimation of parameters that appear in quantum field theory including proper times and accelerations.

C. Sabin, J. Kohlrus, D. E. Bruschi and I. Fuentes, EPJ Quantum Technology 3, 8 (2016) builds on the knowledge that quasiparticles in a Bose-Einstein condensate are sensitive to space-time distortions, and that gravitational waves can induce transformations on the state of phonons that can be observed through quantum state discrimination techniques. This papers shows that this observation method is highly robust to thermal noise and depletion. A bound on the strain sensitivity is derived that shows that the detection of waves in the kHz regime is not significantly affected by temperature in a wide range of parameters that are well within current experimental reach.

It is an object of the current invention to provide an alternative quantum gravimeter that exploits the relativistic effect of phonon creation by gravitational waves.

SUMMARY OF THE INVENTION

According to a first aspect of the invention we provide a method of measuring a gradient of a field (such as a gravity gradient) using one or more trapped BECs, each having a density distribution, and at least two atomic modes within those BECs; wherein the density distributions of the one or more BECs are modified by the field, which leads to a change of the atomic modes' time evolution; wherein the method comprises:

-   -   selecting two atomic modes that are differently affected by the         field; measuring a field-induced difference between the selected         atomic modes; and using the measured field-induced difference to         infer a gradient of the field.

The field may be a gravitational field, and the field induced difference between the two modes may comprise a field-dependent phase. Alternatively, or additionally a number of excitations may be exchanged between the modes or created within one or both of the modes.

The method may comprise using a single trapped BEC having a density distribution that is perturbed by the external field. The two differing modes (e.g. two modes whose field-induced phase differs) may be selected from modes existing in the single BEC.

Each atomic mode may be defined by an angular momentum quantum number and a magnetic quantum number. Measuring the field induced phase difference between two modes may involve selecting two modes with the same angular momentum quantum number, but with differing magnetic quantum numbers. The magnetic quantum numbers of the selected modes may differ by one.

The or each BEC may be trapped in a trapping potential. The trapping potential may comprise an electromagnetic portion and a field-induced portion. The field-induced portion may be due to a gravitational field, which may be caused by a mass distribution in the vicinity of the BEC. The mass distribution may cause a distortion in the shape of the BEC.

The two modes may be at an angle to each other, such that the distortion in the shape of the trapping potential causes the two modes to have differing effective lengths. The two modes may be perpendicular to each other.

The electromagnetic portion of the trapping potential may be a harmonic trapping potential. The gravitational field may add a linear term to the electromagnetic portion of the trapping potential.

A plurality of independent consecutive measurements of the phase difference may be taken. The measurements may then be averaged to produce an averaged phase difference. This averaging may improve the accuracy of the final result.

The one or more atomic modes within the BEC may be prepared in a squeezed state. This may result in a further improvement in the accuracy of the measurement.

According to a second aspect of the invention we provide a quantum gradiometer operable to perform the method of the first aspect of the invention.

According to a third aspect of the invention we provide a system for measuring a field gradient; the system comprising:

-   -   a trapped BEC having a density distribution, wherein the trapped         BEC comprises a plurality of atomic modes, and wherein the         density distribution of the trapped BEC is modified by the         field;     -   a measurement system for measuring a field-induced difference         between two selected atomic modes within the trapped BEC; and     -   a processor for inferring a field gradient from the measured         field-induced difference.

The field induced difference may be a phase difference between the modes, a number of excitations exchanged between the modes, or a number of new excitations created within one or both of the modes.

The system may be operable to perform the method of the first aspect of the invention, including optional features.

According to a fourth aspect of the invention we provide a method of measuring an acceleration due to a gravitational field using one or more BECs,

-   -   wherein the one or more BECs include:         -   a first mode propagating in/parallel to a first direction             (e.g. between a first location (which may be a boundary) r             _(b) and a second location (which may be a boundary) r             _(t)), wherein the first direction may be parallel to a             direction of the acceleration due to the gravitational             field; and         -   a second mode propagating in/parallel to a second direction,             different to the first direction, between a third location             (which may be a boundary) z_(l) and a fourth location (which             may be a boundary) z_(r), wherein the second direction may             be perpendicular to the first direction and to the direction             of the acceleration due to the gravitational field;     -   wherein the method comprises:         -   measuring a phase difference between the first mode and the             second mode; and         -   using the measured phase difference to measure the             acceleration due to the gravitational field.

Such a scheme is advantageous as is makes use of collective oscillations (i.e. modes) within the at least one BEC, rather than relying on atom interference. It can thus be made smaller than prior atom interferometers.

The first direction is, in most realistic scenarios, a vertical direction, and the second direction is thus a horizontal direction. The frequency in the vertical mode depends on the Schwarzschild radius, as the acceleration due to gravity causes a density perturbation in the BEC in the vertical direction, which has an effect on the frequency of the vertical mode. A phase difference between the modes thus can be used to provide information about the Schwarzschild radius (and therefore gravity).

The phase difference between the first frequency mode and the second frequency mode may be used to estimate the local Schwarzschild radius r_(s).

The local gravitational field may be estimated from the local Schwarzschild radius r_(s) using

g=r _(s) c ²/2 r _(b) ²

-   -   r _(b)

The one or more BECs may have a uniform density when not under the influence of an external potential.

The method may further comprise condensing the at least one BEC in a trapping potential which is tuned to partially cancel the acceleration caused by the gravitational field, such that a density perturbation δ_(ρ) is created in the BEC.

The trapping potential may be given by

$\mspace{76mu}{{\Delta\; V} \approx {\frac{m^{2}c^{2}}{\text{?}^{2}}\frac{r_{S}x^{\text{?}}}{{\overset{\_}{r}}_{b}^{2}}}}$ $\mspace{76mu}{\overset{\_}{r}}_{b}$ ?indicates text missing or illegible when filed

and the corresponding relative density perturbation δ_(ρ) is given by

$\delta_{\rho} = {\frac{1}{2}\frac{c^{2}}{c_{s\; 0}^{2}}\left( {\frac{r_{S}}{{\overset{\_}{r}}_{b}^{2}} - \frac{r_{S,0}}{{\overset{\_}{r}}_{b,0}^{2}}} \right)\left( {x^{f} - \frac{x_{t}^{f}}{2}} \right)}$ ${\overset{\_}{r}}_{b}$ ${\overset{\_}{r}}_{t}$

Alternatively, the at least one BEC may comprise two BECs, and the first (vertical) mode may propagate in the first BEC, whilst the second (horizontal) mode may propagate in the second BEC.

The method may further comprise the steps of:

-   -   condensing a second BEC in a second atom trap having a second         length L_(p) between the third and fourth locations (which may         be boundaries) z_(l) and z_(r), the second length being equal to         the first length,     -   wherein the first mode propagates in the first BEC and the         second mode propagates in the second BEC.

The or each BEC may be confined in a box potential. For instance, the first BEC may be confined in a first box potential and the second BEC may be confined in a second box potential. The first box potential may be effectively one-dimensional, with the mode being used in the analysis propagating along that one dimension. Similarly, the second box potential may be effectively one-dimensional, with the mode being used in the analysis propagating along that one dimension.

The length L_(p) of each of the first and second BECs may be held constant by a respective rigid rod, and deformation of the rigid rod due to gravity may be taken into account for the first BEC. Such a deformation may be neglected for the second BEC.

The deformation of the rigid rod may affect the frequencies of the frequency modes through

ω_(v,k) ^(def)=ω_(h,k)(1+Θ)

where Θ is a frequency shift.

The first mode may be a squeezed state. The method may further comprise creating a squeezed state for a first mode k of the first BEC.

The second mode may be a coherent state. The method may further comprise preparing a coherent state a second mode of the second BEC. The second mode of the second BEC may have the same mode number k as the squeezed mode of the first BEC.

The method may further comprise entangling the first and second modes.

A phase shift may be introduced between the first and second modes. The phase shift may be π/2.

The or each BEC may be modelled as a quantised field.

The first mode may be a phonon mode. The second mode may also be a phonon mode.

The spacetime of the or each BEC may be described using Schwarzschild geometry in order to calculate the acceleration due to gravity.

The phase shift may be used to determine the Schwarzschild radius using

$\begin{matrix} {\mspace{76mu}{\omega_{h,k}^{def} = {\frac{\text{?}\;\pi\; k}{L_{p}}{f\left( r_{b} \right)}^{1/2}\mspace{14mu}{and}}}} & (61) \\ {\mspace{76mu}{\omega_{v,k}^{def} = {{{\omega_{h,k}\left( {1 + \Theta} \right)}.\text{?}}\text{indicates text missing or illegible when filed}}}} & (62) \end{matrix}$

where v and Θ are defined for the light cavity and the BEC respectively as

$\begin{matrix} {{Light}\left\{ {\begin{matrix} {v = c} \\ {\Theta = \Lambda} \end{matrix},{{BEC}\left\{ {\begin{matrix} {v = c_{s\; 0}} \\ {\Theta = {{\frac{3}{2}\Lambda} - \frac{\delta_{\rho}^{r\; 2}L_{p}^{\; 2}}{64}}} \end{matrix}.} \right.}} \right.} & (63) \end{matrix}$

According to a fifth aspect of the invention we provide a quantum gravimeter comprising:

-   -   r _(b) r _(t)     -   a second BEC in a second atom trap having a second length L_(p)         between the third and fourth locations z_(l) and z_(r),     -   r _(b) r _(t)     -   a processor operable to use the phase difference between the         first mode and the second mode to measure the acceleration due         to the gravitational field.

According to a sixth aspect of the invention we provide a quantum gravimeter comprising:

-   -   r _(b) r _(t)         is perpendicular to the direction of the acceleration due to the         gravitational field and to the first longitudinal axis;     -   r _(b)     -   r _(t)     -   a processor operable to use the phase difference between the         first mode and the second mode to measure the acceleration due         to the gravitational field.

It will be appreciated that the features described above with respect to the method of the fourth aspect of the invention may be implemented in a quantum gravimeter according to the fifth aspect of the invention or the sixth aspect of the invention.

According to a seventh aspect of the invention we provide a quantum gravimeter/accelerometer operable to perform the method of the fourth aspect of the invention.

The technology described herein can be adapted to measure any acceleration, not just the acceleration due to gravity, and therefore according to an eighth aspect of the invention we provide a method of measuring an acceleration (e.g. that experienced when at rest in a linear gravitational potential) using one or more trapped BECs and at least two atomic field modes within those BECs; wherein the density distributions of the one or more BECs are modified by the acceleration, which leads to a change of the field modes' time evolution; wherein the method comprises:

-   -   selecting two modes that are differently affected by the         acceleration; and     -   measuring an acceleration-induced difference; and     -   using the measured acceleration induced difference to infer the         acceleration.

The acceleration-induced difference may comprise one or more of:

-   -   (i) an acceleration dependent phase;     -   (ii) a number of excitations exchanged between the modes;     -   (iii) a number of new excitations created within one of the         modes, or both of the modes.

The two modes may be at an angle to each other, such that the acceleration causes the two modes to have differing effective lengths. The two modes may be perpendicular to each other, for example, the two modes may comprise a mode that is parallel to a direction of the acceleration and a mode that is perpendicular to a direction of that acceleration.

The acceleration may be an acceleration due to gravity, and thus the method of the eighth aspect of the invention may be combined with any of the method steps of the fourth aspect of the invention if required.

The above methods may also be used to measure/analyse the field gradient (e.g. a gravity gradient) rather than the acceleration due to a field. Density variations in one or more BECs caused by a field can be used to infer a local gradient of that field by measuring a phase difference between two differently oriented modes within the one or more BECs. The techniques of the 4^(th)-8^(th) aspects of the invention may thus be combined with the techniques of the 1^(st)-3^(rd) aspects of the invention.

BRIEF DESCRIPTION OF THE DRAWINGS

Aspects of the invention will now be described in more detail, by way of example only, with reference to the accompanying drawings, in which

FIG. 1 schematically illustrates a quantum gravimeter;

FIG. 2 schematically shows a) an effective uniform Bose-Einstein condensate (BEC) confined within a 1-dimensional box potential, and b) as a comparison for the BEC case, an effectively 1-dimensional light field confined using mirrors;

FIG. 3 schematically illustrates effects of acceleration on a rod used to suspend the device creating the trapping potential for the BEC;

FIG. 4 schematically illustrates the principle of an interferometric measurement set up;

FIG. 5 sets out steps of a method of measuring acceleration; and

FIG. 6 sets out a schematic proposal for a detector

FIGS. 7a and 7b show the effects of a gravitational mass distribution on the BEC

FIG. 8 illustrates a scheme for measuring density of a BEC;

FIG. 9 schematically illustrates a quantum gradiometer; and

FIG. 10 sets out steps of a method of measuring a field gradient.

DETAILED DESCRIPTION

We have developed a measurement scheme that changes the paradigm on how field effects, particularly gravitational effects, are measured using quantum systems. Our scheme overcomes limitations due to time of flight constraints by using the effects of a field, such as gravity, on collective excitations (called phonons) of static and oscillating Bose-Einstein condensates (BECs). Two exemplary measurement schemes are described below. First, a system and method for measuring an acceleration due to a field (e.g. gravity) is summarised. Second, a system and method for measuring a field gradient (e.g. a gravity gradient) is summarised. Both methods are then described mathematically in more detail in the sections that follow.

Gravimeter

An example system for measuring an acceleration due to a field is schematically illustrated in FIG. 1 as gravimeter 10. The gravimeter 10 may be used to very accurately measure an acceleration due to a gravitational field (proper acceleration in this case matches gravitational acceleration).

The gravimeter comprises at least one, and preferably first and second atom traps 12, 14, which are perpendicular to one another. The first atom trap 12 is aligned with a direction of an external acceleration to be measured, in this case acceleration due to gravity, indicated in FIG. 1 by the arrow g. The first atom trap is thus aligned along a direction that is radial to a surface of a body causing the gravitational acceleration which is to be measured (e.g. the Earth), and so for most purposes the first atom trap is substantially vertical. The second atom trap 14 is perpendicular to the first atom trap 12, and hence perpendicular to the acceleration being measured. For most purposes the second atom trap is substantially horizontal.

Each of the atom traps 12, 14 contains, when in use, a Bose Einstein condensate (BEC) confined within a trapping potential. A BEC is a state of matter comprising bosons, in which a large fraction of the bosons occupy the lowest quantum state. A BEC may be formed by cooling a dilute gas of bosons below a critical temperature, very close to absolute zero (−273.15° C.). Methods of forming BECs are well known in the art, and so will not be discussed in further detail herein.

-   -   r _(b) r _(t)

Measuring equipment 34 is operable in use to measure a phase difference between first frequency of a first mode of the first BEC 26 and a second frequency of a second mode of the second BEC 30. A processing device 36 is operable to infer information about the local gravitational field from the frequency difference, and in particular from the difference in the modification of the frequencies of the modes due to the difference in the modification of the BEC density induced by acceleration/gravity, as explained in more detail below.

In particular, the gravimeter may be used to measure the Schwarzschild radius (ie. the product 2GM/c² where G is the gravitational constant and M the mass) due to a gravitational field. The gravimeter relies on using the effects produced by the spacetime of the Earth (or by gravity) on the quantum states of a phononic quantum field and the bulk (i.e. the atoms in the ground state, also referred to as “atomic density”) of a BEC. These effects include a perturbation in the atomic density of the vertically oriented BEC that is due to the gravitational field.

In previous quantum gravimetry methods atom-atom interactions are considered as a systematic effect. In contrast, in this method we use two-body collisions that give rise to entanglement between atoms and matter wave vibrational modes (called “phonons”, which correspond to excitations also referred to as “density fluctuations”). In the absence of collisions there are no phonons. The reason why we use two-body interactions is because the dynamical evolution of phononic modes that they produce is affected by gravity.

The gravitational field affects the frequencies of these modes through the quantised field equation

$\begin{matrix} {{{{\hat{\psi}}_{h}\left( {x^{0},\overset{\sim}{z}} \right)} = {\sum\limits_{k}\left\{ {{{\hat{a}}_{h,k}{\phi_{h,k}\left( {x^{0},\overset{\sim}{z}} \right)}} + {{\hat{a}}_{h,k}^{\dagger}{\phi_{h,k}^{*}\left( {x^{0},\overset{\sim}{z}} \right)}}} \right\}}},} & (24) \end{matrix}$

The solutions to this equation are plane waves ϕ_(h, k) (k indicates the mode number) that have frequencies which depend on the height in the gravitational field x, the Schwarzschild radius r_(S)=2GM/c² (where c is the speed of light), the length of the BEC Lp and the effective speed of sound c_(s) (ie. the atomic mass, density and the strength of the interactions (scattering length). Information about the local Schwarzschild radius r_(S) (and thus the local gravitational field) can therefore be obtained by analysing the frequencies of phonon modes within a BEC.

The Schwarzschild radius is defined as:

r _(S)=2GM/c ²  (A)

The local gravitational field g can be obtained using

g=GM/r ²  (B)

Replacing GM in equation (B) gives

g=r _(s) c ²/2r ²  (C)

where r is the radial height coordinate within the gravitational field. The local gravitational field can thus deduced from the local Schwarzschild radius using equation (C).

Further details explaining why there is field information encoded in the difference between the frequency modes, and how that information may be extracted, is provided in the detailed discussion below (headed DETAILED DISCUSSION: GRAVIMETER). First, however, we give further information on an example experimental set-up.

An example of a one-dimensional trap 16 suitable for use in the gravimeter of FIG. 1 is the “box trap” shown in FIG. 2. The box trap is formed using one hollow tube laser beam 18 and two sheet laser beams 20, 22, which together create a repulsive potential 16 for the atoms within the trap. The atomic cloud is confined by the three beams to the dark cylindrical trapping region 24, which can be considered effectively one-dimensional for our purposes in that the trapping region 24 is very small in two dimensions and much longer in a third dimension, which corresponds to the length Lp of the trapping region 24. Such a set-up allows us to assume that the BEC is only excited for momenta in the longitudinal dimension (i.e. the dimension having length Lp). It will be appreciated that there are multiple ways to make such a potential, and this is an example only. Other types of potential (e.g. a harmonic potential) could also be used, if required.

An atom trap of the type discussed above can be created using 532 nm laser beams. With a total laser power of P₀≈700 mW it is possible to achieve a trap depth of V₀≈k_(B)×2 μK (where k is the Boltzman constant). Further details of the construction of an appropriate atom trap can be found in A. L. Gaunt, T. F. Schmidutz, I. Gotlibovych, R. P. Smith, and Z. Hadzibabic, Phys. Rev. Lett. 110, 200406 (2013).

A first BEC 26 may be condensed within the first atom trap 12 using standard cooling techniques. Similarly, a second BEC 30 may be condensed within the second atom trap 14 using standard cooling techniques.

Accuracy may be improved by, for the vertical BEC, creating a squeezed state for a mode k (the method works with other states but is improved by using squeezed states). Squeezing phonons can be done by several methods that are well known in the art (such as changing the length of the trap, adding an electromagnetic potential, etc).

For the horizontal BEC, a mode within the same mode number k can be prepared in a coherent state (although a squeezed or other state can also be used if preferred). Note that the horizontal mode can have a different mode number, but it is more convenient mathematically to use the same mode number. Indeed, the horizontal “mode” could be the bulk of the BEC rather than a phonon mode, which may result in greater sensitivity (since the number of atoms can be much higher in the bulk, and the bulk is usually in a coherent state).

As explained above, the gravimeter is oriented with the first longitudinal axis aligned with the direction of the acceleration due to the gravitational field (i.e. vertical/radial), and the second longitudinal axis perpendicular to the first (i.e. horizontal, or alternatively tangential). Information about the gravitational field is encoded within the vertical mode, but not in the horizontal mode (as will be explained below)

The horizontal BEC is located adjacent an end of the vertical BEC. In the example shown the horizontal BEC 30 is located adjacent a lower end of the vertical BEC 26. In particular, the longitudinal axis 32 of the horizontal BEC is located at the same vertical coordinate as the lower boundary r _(b) of the vertical BEC. Alternatively, the horizontal BEC may be located at the upper boundary r _(t) of the vertical BEC. As explained in more detail below, the accuracy of the gravimeter is increased when the horizontal BEC is located at an end of the vertical BEC; however other locations are possible if required, with the exception that the horizontal BEC may not be located halfway along the vertical BEC (i.e. at r _(b)+L_(p)/2)—further details as why this is the case are provided below.

Each BEC 26, 30 is confined within an atom trap which confines the atoms within the trap in a trapping potential V₀ created by the box trap arrangement of lasers. The gravitational field of the Earth affects the atomic density in the vertical box potential V₀, such that it will not remain uniform unless the effect due to the Earth's gravitational field is cancelled out. We can assume that the gravitational field has no effect on the horizontal box potential V₀ because the potential is effectively one-dimensional and perpendicular to the direction of the gravitational field, meaning the horizontal potential experiences no force due to the gravitational field in its longitudinal dimension.

A linear electromagnetic potential ΔV is applied to the vertical trap 12 (for example by applying a magnetic field of the right strength) in order to cancel out the effects of the external gravitational field. In particular the linear potential ΔV is applied to the vertical BEC 26 that cancels the effect of the gravitational field on the bulk of the atoms within the BEC. The field is typically cancelled up to second order in r_(s) √{square root over (r_(b))} and x√{square root over (r_(b))} (where x is the radial (i.e. vertical) coordinate and r _(b) is the position of the lower boundary in isotropic Cartesian coordinates—discussed in more detail below), and up to the first order in their product. Therefore the density will be uniform to that order. In practice, an experimentalist may tune the gradient of the linear potential dΔV/dx such that the density is homogeneous.

Alternatively, we have found that dΔV/dx can also be fixed to a value that differs slightly from the value necessary to cancel the full effect of the gravitational field. This difference leads to a density perturbation in the first BEC 26. We can use this perturbation to increase the accuracy of our estimate of the Earth's spacetime parameters, as explained in more detail below.

The accuracy of the estimate may also be increased by taking into account the effects of length distortion in the vertical direction due to the gravitational field. Again, this effect may be disregarded in the horizontal direction.

Many techniques exist for the measurement of the atomic mode frequencies, for instance counting the phonon difference in each mode (vertical and horizontal), dual parameter estimation, parity measurement or using Bayesian protocol.

A variant of the method is that instead of using two separate one dimensional BECs, one could instead use a single BEC in a three-dimensional box potential and use two modes in that single BEC, one horizontal mode in a coherent state and one vertical mode in a squeezed state

It will be appreciated that the two modes do not need to be vertical and horizontal, nor do they need to be perpendicular to one another. The arrangement shown in FIG. 1 is useful in that the mathematics may be simplified (as described more below) but other modes may be used if required. Any two modes that are differently affected by the field being measured may be selected if required.

A third variant of the method, which can also be implemented using either two separate BECs or a single BEC in a three dimensional box involves transferring excitations between the bulk and the two selected phononic modes (e.g. vertical and horizontal) using a beamsplitter gate and extra phase shift.

Many techniques exist for the measurement of the output of 30, for instance counting the phonon difference in each mode (vertical and horizontal), dual parameter estimation, parity measurement or using Bayesian protocol.

A fourth variant of the method, which can be implemented using either two separate BECs or a single BEC in a three dimensional box, is a pumped SU(1,1) scheme which involves transferring excitations between the bulk and the phononic modes using for example tritter gates and parametric amplification. Other gates mixing excitations between the bulk and phononic modes can be also used

Contrary with atom interferometry where interactions are avoided, the gravimeter described above takes into account two-body collisions.

In the example implementation discussed above the BEC measures 100 μm, therefore, our design can be implemented with state of the art table top technologies (including laser system and vacuum chamber) but soon this technology will be miniaturised into chips. The phononic field can be cooled down to 0.5 nK making the device highly sensitive.

Current quantum gravimeters are based on non-relativistic quantum mechanics and cannot be employed to measure gravity beyond the Newtonian regime because quantum theory and general relativity are incompatible theories. This theoretical framework enables the description of both quantum theory and general relativity in a consistent way. The detection scheme relies on the deformation of the BEC (a relativistic effect) produced by the gravitational field of the Earth on a Bose-Einstein condensate consisting of atoms cooled below the nano-Kelvin regime. Vertical deformation produces an observable effect on the phononic modes of the BEC as compared to longitudinal modes. These differences between the vertical and longitudinal modes are observable and can be used to estimate the local gravitational field.

Note that BECs are also used in atom interferometry for high precision experiments. However in atom interferometry, the wave function of single atoms is brought into interference with itself; no collective phenomenon is used, while phonons in BECs are the quasi-particles of collective oscillations of atoms.

Such a gravimeter can be used to measure the proper acceleration and/or the Schwarzschild radius of the Earth. As discussed in significantly more detail below, the basis of the measurement is the difference of the fundamental frequency of phonon modes in differently oriented BECs (or in differently oriented modes within a single BEC) due to the deformation of the length of the trapping potential of the BEC (or the length of the cavity). The relative error bound for the estimation of the Schwarzschild radius on the length scale of 200 μm is 10⁻⁸ in one run (6 seconds) and 10⁻¹¹ if one integrates over a year. This compares very well with commercial atomic quantum gravimeters that reach relative errors of 10⁻⁹ in 15 days with a device that measure 100×50×70 cm³. The gravimeter discussed herein has the advantage that the device can be miniaturised without losing precision and can be employed to measure field strengths beyond the Newtonian regime.

Gradiometer

Another application for BECs is in gradiometry, that is, to measure differential forces due to a field gradient, such as a gravity gradient.

FIG. 6 shows a basic scheme for a ground state interferometer. A BEC is split into two clouds and then recombines to interfere. The interference measures any field gradient due to an object and constrains deviations from general relativity (e.g. fifth force models). Such a detector has a long integration time: atoms can be held near the object for a longer time than with a ballistic scheme.

However, the measurement sensitivity can be improved using collective quantum BEC excitations, as demonstrated in the present invention. An initially squeezed phonon state has better phase estimation and error scaling than any classical state.

FIG. 7a shows a BEC 38 trapped in a potential, which in this case is a spherically symmetrical, harmonic trapping potential 42. Such a potential may be generated by light and/or magnetic fields in a known manner. A BEC 38 in such a trapping potential has a known initial spatial density distribution, which in this case is a density distribution that varies parabolically with distance from the centre of the trap. The BEC can thus be considered to be a sphere having a substantially constant radius

o. Atomic field modes can exist with such a trapped BEC. Atomic field modes (also referred to as phonons) are collective oscillations of the atoms in a BEC 38 with well-defined (quantised) energy levels.

If an arbitrary mass distribution 46 is placed in the vicinity of the BEC 38, the BEC 38 will deform due to the gravity exerted by that mass. This deformation changes the properties of the BEC, in particular the atomic mode structure of the phonon field.

Typically, using the present scheme one can measure a 20 mg mass at the distance of 1 mm with a good signal to noise ratio. For every three orders of mass increase, one gets one order of magnitude increase in distance (e.g. 20 kg at a distance of 1 m, and 20 tonnes at a distance of 10 m).

A gradient of the field causing the deformation (i.e. a gradient of the gravitational field) may be determined by measuring a phase difference between two atomic field modes that are differently affected by the field. Such modes may be, for example a first mode that is at a first orientation (e.g. parallel) with respect to a selected axis along which a field due to the mass distribution acts, and a second mode that is at a second, different, orientation (e.g. perpendicular) to that axis. As the BEC is extremely cold there are only a few thermal excitations of the phonon field, which represent the noise background on which we measure.

FIG. 7b shows the BEC 38 after it has been deformed due to a mass distribution 46 in the vicinity of the BEC 38, in this case a spherically symmetric, homogenous mass distribution. The mass distribution causes a force to act on the BEC along a known axis (shown in this example as the z-axis). The presence of the mass distribution 46 thus adds a gravitational (Newtonian) potential to the original trapping potential 42, resulting in a new trapping potential 44 with its centre shifted by amount z_(g) along the z-axis, as shown in FIG. 7b . As a result, the centre of mass of the BEC 38 is also shifted by amount z_(g) along the z-axis. Furthermore, the trapping potential 44 profile is also broadened along the z-axis, when compared to the profile of the original trapping potential 42. The overall effect is that the BEC 38 is elongated (or deformed) along the z-axis, and assumes a spheroidal shape having a second, deformed, density distribution (which is in this example also a parabolic distribution). The radius of the trapped BEC thus can no longer be considered substantially constant, but instead varies between a maximum radius

/λ (which is along the direction parallel to the field due the mass distribution, i.e. the z-axis) and a minimum radius

perpendicular to the direction of that field).

The spheroidal shape and the parabolic density distribution of the BEC 38 give rise to a well-defined set of atomic field modes within the BEC that are different to the field modes in an unperturbed BEC. Those perturbed atomic field modes, and in particular the angular frequencies/wavelengths of those modes, encode information about the gravity gradient.

One approach to determine gradient information about the field that is perturbing the BEC from a perturbed BEC is to select any two arbitrary modes of the perturbed BEC 38. The two modes can be allowed to evolve for a fixed time interval, and then the phase difference can be measured experimentally to reveal the gravity gradient. Therefore, one of the modes effectively serves as a reference for the other. As discussed in more detail below, such a system works effectively if the two selected modes have the same total angular momentum quantum number, but have magnetic quantum numbers that differ (e.g. by one).

In the specific scheme illustrated, of a BEC in a harmonic trap, information about the field gradient can thus be determined by considering phonons with two specific angular momentum quantum numbers. Any two phonons modes with different angular momentum quantum numbers may be selected. In one example, we select modes with angular momentum projected in the z-direction l and l−1, where l is the absolute (total) value of the angular momentum of both modes. For these modes, the frequencies are known when the BEC is deformed (ellipsoidal) in the z-direction. These frequencies in the gravitational field of a sphere are:

ω_(l,l)≈√{square root over (l)}(ω₀ +GM/2R ³ω₀).

and

ω_(l,l−1)≈√{square root over (l)}[ω₀+(l−3)GM/2lR ³ω₀].

∈_(grad):=2MG/R ³

ϕ=ωt

Then, the principle limit on the precision of a measurement of the gravity gradient through this phase is given by

$\begin{matrix} {{\delta_{e_{grad}} = {\frac{2\omega_{0}^{2}}{{a(l)}\epsilon_{grad}}\frac{1}{\sqrt{2{{Nn}\left( {n + 1} \right)}}}}},} & (2) \end{matrix}$

A specific way to access the phase of the phonon modes could be let these modes interfere. This is described in equation (3) giving the frequency difference between the two modes:

$\begin{matrix} {{\Delta\omega}_{l} = {{{\omega_{l,l} - \omega_{l,{l - 1}}} \approx {\frac{3_{\epsilon_{grad}}}{4\sqrt{l}\omega_{0}}.{\Delta\phi}_{l}}} = {{\Delta\omega}_{l}{t.}}}} & (3) \end{matrix}$

can obtain an absolute sensitivity of 10⁻⁸ s⁻² for gravity gradient measurements using this scheme.

However, it can be appreciated that the present invention is not limited to measuring only the phase difference of the two phonon modes whose frequencies are dictated above—other modes can be selected. It can also be appreciated that not all phonons modes may have an analytical representation. In such cases, numerical methods may be employed to order to obtain the relevant parameters, such as the frequencies of the modes.

There exist many methods of extracting the phase difference. For example, in order to experimentally measure the phase difference, one can either try to measure the phase difference directly by allowing the modes to interfere, as mentioned above, or to infer the phase difference of the modes separately by measuring, for example, the density of the BEC 38 and then calculating the phase difference from this.

One method of calculating the phase difference between phonon modes by measuring the density of a BEC comprising the modes has been termed ‘heterodyne detection’ by some authors, and is illustrated in FIG. 8. After the modes are excited in the BEC 38, the trapping potential is switched off and the BEC starts to expand and fall freely under gravity. During the expansion the energy contained in the phonons is transformed into the kinetic energy of atoms. These particles may interfere with the atoms in the ground state within the BEC 38 to produce periodic density modulations which can be measured using a detector 48. In the present invention, the periodic density modulations can result from interference of the atoms that have derived their kinetic energy from two different modes of the BEC. Thus, these density modulations contain the information about the phase difference of the phonons. Information about the local gravitational field gradient can be inferred from the phase difference, as discussed above.

It is also possible to control the phonon-phonon interactions in the BEC with Feshbach resonances via manipulation of the BEC's shape, enabling multi-mode operations that can be used for the implementation of Mach-Zehnder interferometer (MZI) schemes. The controlled coupling of the phonons is equivalent to a beam splitter in the MZI, and thus, can be used to infer the phase or phase difference of the phonons directly.

Typically, the measurements made using the above methods are carried out multiple times. The measurements are then averaged to produce an average phase difference, to improve accuracy of the final result.

Further improvement in accuracy is obtained if phonons are created in a squeezed state. Squeezing could be achieved by periodically changing the shape of the BEC (for example by modulating the trapping potential) or the atom-atom interaction strength of the BEC.

FIG. 9 shows an exemplary gradiometer which implements the methods discussed above. A BEC 40 is trapped in a trapping potential and is deformed by an external field acting in the direction F. Measuring equipment 48 is operable in use to measure a phase difference between a first mode of a perturbed BEC 38 and a second differently oriented mode of the same BEC. A processing device 50 is operable to infer information about the gradient of the field perturbing the BEC from the measured phase difference.

As an alternative to the schemes discussed above, two BECs could be used (similar to the gravimeter scheme discussed above), where one of the BECs is perturbed by the field whilst the other BEC is not perturbed, or is differently perturbed. A method and system for extracting the gravity gradient information would then involve measuring the phase differences between two modes, one ‘reference’ mode in the unperturbed BEC and another ‘probe’ mode in the perturbed BEC, both of which are allowed to evolve for a fixed time interval. The reference mode may be a selected mode k of the unperturbed BEC, and the probe may be the same mode k (i.e. with the same angular momentum quantum number and magnetic quantum number), but of the perturbed BEC. As explained above, the phase difference could be extracted by allowing the modes to interfere.

Detailed Discussion: Gravimeter Section I. Summary of Detailed Analysis

The theoretical framework underpinning the gravimeter described above is discussed in the pages below. Besides the phonon field in a BEC, we consider a light field in a 1-dimensional optical cavity to compare with. In particular, we study the frequency spectrum of the quantum fields when they are kept at a fixed height above the ground in a laboratory on the surface of the earth. This situation corresponds to a proper acceleration that matches the gravitational acceleration. We consider two orientations of the BEC/cavity, horizontal and vertical, that lead to different effects on the frequency spectrum. The time evolution of any mode of the light field/phonon field can be used as a clock. We derive bounds for the precision of these extended clocks and discuss the effect of time dilation and proper acceleration. The influence of time dilation on the evolution of the extended clocks can be used for the estimation of gravitational red shift and the Schwarzschild radius of the earth by comparing the readings of two clocks at different height. We give limits for the precision of such measurements. The proper acceleration of the extended clocks matches the gravitational acceleration. Therefore, the different effect of proper acceleration on the clocks in horizontal and vertical orientation can be used to infer the Schwarzschild radius of the Earth locally. We derive the principle limitation for such a measurement. We also present an experimental setup that can be used to approach the principle bound for the local estimation of the Schwarzschild radius: modes of the quantum field in differently oriented BECs/cavities are brought into interference in a Mach-Zehnder interferometer (MZI) type setup.

The description below is organised as follows: in Section II we present the physical set-up; a hovering 1-dimensional BEC/cavity in the gravitational field of the Earth. In Section III, we introduce our model of a 1-dimensional BEC in the gravitational field of the Earth and we introduce the phonon field which serves as the counterpart to the light field in an optical cavity, as well as the deformation of the density due to gravity. The quantisation of both the phonon field in a BEC and the light field in a cavity, is presented in Section IV and the frequency spectrum is derived for horizontal and vertical orientation. In Section V, we derive the effect of gravity on the length of the BEC/cavity in the vertical orientation. This leads to an additional change of the frequency spectrum that we compare with the relativistic effect. Additionally, we consider the cavity as an extended clock and we apply methods of quantum metrology to provide bounds on the precision for the estimation of proper time. Furthermore, we give the precision for the estimation of red shift and the Schwarzschild radius by independent measurements at different heights. In Section VI, we introduce local estimations of gravitational accelerations, i.e. the Schwarzschild radius, with two differently oriented BECs/cavities. Finally, in Sec. VII we propose a Mach-Zehnder interferometer-like setup, which can in principle approach the precision bound. A summary, conclusion and an outlook on possible applications and future work are given in Section VIII.

The conventions used herein are the following: the metric signature is (−, +, +, +) in 3+1-dimensions and (−, +) in 1+1 dimensions. We use the symbol Tp for transposition when convenient to distinguish it from time.

Section II. The Hovering I-Dimensional Cavity

Effects of earth rotation and deviations from spherical symmetry can be safely ignored for our purposes. Therefore, we describe the spacetime surrounding the earth using the Schwarzschild geometry, which is given as

$\begin{matrix} {{g - {{diag}\left( {{- {f(r)}},\frac{1}{f(r)},r^{2},{r^{2}\sin^{2}\theta}} \right)}},} & (1) \end{matrix}$

in spherical Schwarzschild coordinates (x^(o), r, υ, φ). Here, x^(o)=ct, f(r)=1−r_(s)/r, the Schwarzschild radius is given as r_(s):=2GM/c², the parameter M is the mass of the planet and G is the gravitational constant. For our purposes it is useful to consider the transformation of the Schwarzschild metric in Eq. (1) to the isotropic Cartesian coordinates (x^(o), x, y, z). With r=(1+r_(S)/4r)² r where r=(x²+y²+z²)^(1/2), we arrive at the Schwarzschild metric in isotropic spherical coordinates. After a the usual transformations from spherical to Cartesian coordinates, we find

$\begin{matrix} {\mspace{76mu}{g_{\text{?}} = {\left( {1 + \frac{r_{S}}{4\overset{\_}{r}}} \right)^{4}{{{diag}\left( {{- \frac{\left( {1 - \frac{r_{S}}{4\overset{\_}{r}}} \right)^{2}}{\left( {1 + \frac{r_{S}}{4\overset{\_}{r}}} \right)^{6}}},1,1,1} \right)}.\text{?}}\text{indicates text missing or illegible when filed}}}} & (2) \end{matrix}$

Note that the coordinate transformation from the Schwarzschild coordinates to the isotropic coordinates is time independent. Therefore, any point at a fixed coordinate position in Schwarzschild coordinates stays at a fixed coordinate position in isotropic coordinates. Hence, there is no additional proper acceleration introduced by fixing the coordinate position in the new coordinate system.

$\mspace{70mu}{r_{S}/\overset{\_}{r}}$ $\mspace{70mu}{r_{S}/\overset{\_}{r}}$ $\begin{matrix} {{{g_{\mu\nu} = {\eta_{\mu\nu} + {\frac{r_{S}}{\overset{\_}{r}}{{diag}\left( {1,1,1,1} \right)}} + {\frac{r_{S}^{2}}{2{\overset{\_}{r}}^{2}}{{diag}\left( {{- 1},{3/4},{3/4},{3/4}} \right)}}}},\mspace{70mu}\frac{r_{S}^{\text{?}2}}{{\overset{\_}{r}}^{2}}}\mspace{70mu}\frac{r_{S}^{\text{?}2}}{{\overset{\_}{r}}^{2}}{\text{?}\text{indicates text missing or illegible when filed}}} & (3) \end{matrix}$

The system of interest is a massless scalar field ψ trapped inside a rectangular cavity. We assume that the cavity is very small in two dimensions and longer in the third one. Such a set-up can describe both, the electromagnetic field confined within an optical cavity or the phonon field inside a trapped BEC. In the latter case, the bulk of the BEC acts as a cavity for the phonons. In the following, we introduce the formalism in full generality and we then specialise to the phonon excitations in a trapped BEC.

We assume that the field is only excited for momenta in the elongated dimension and it is found, with good approximation, in the ground state with respect to the transversal dimensions. This assumption can be understood by considering the energy spectrum of the transverse dimensions; the energy difference between two states is inversely proportional to the length of the relevant dimension. The smaller the confinement, the larger the energy difference. Therefore, we can consider the field as effectively 1-dimensional.

-   -   r _(b) r _(t).     -   r _(b)

The edges of the cavity, as well as every other component of the cavity, do not follow geodesics and therefore must have an outward pointing four-acceleration. We anticipate that this will result in a measurable difference between the time evolution of the fields in the vertical and the horizontal orientation.

$\mspace{76mu}{L_{p} = {\int_{a}^{b}\ {\overset{\text{?}}{d}\;\sigma\overset{\text{?}}{\sqrt{g_{\mu\nu}{s^{\mu}(\sigma)}{s^{\nu}(\sigma)}}}}}}$ $\mspace{76mu}{\overset{\_}{r}}_{b}$ $\mspace{76mu}{{{\overset{\_}{r}}_{t}.\text{?}}\text{indicates text missing or illegible when filed}}$

For the horizontal orientation, we start with the coordinates (x^(o), z), and we write r(z):=(r _(b) ²+z²)^(1/2). We consider distances of at most the length of the trap/cavity L_(p) which is of the order between centimeters and micrometers. Therefore L_(p)/r _(b) is approximately of the same order as r_(S)/r _(b). Expanding the metric (3) in second order in s/r _(b) and r_(S)/r _(b) and in first order in their product, we find

g _(h,μv) =f _(RS)( r _(b))diag(−1,Σ( r _(b))),  (4)

where we defined the red shift function:

$\begin{matrix} {{f_{BS}\left( \overset{\_}{r} \right)}\text{:=}\mspace{14mu}{\left( {1 - \frac{r_{S}}{\overset{\_}{r}} + \frac{r_{S}^{2}}{2{\overset{\_}{r}}^{2}}} \right).}} & (5) \end{matrix}$

and the spatial scale function:

$\begin{matrix} {\mspace{79mu}{{\sum\left( \overset{\_}{r} \right)}:={{\left( {1 + \frac{2r_{S}}{\overset{\_}{r}} + \frac{15r_{S}^{2}}{8{\overset{\_}{r}}^{2}}} \right).\mspace{79mu}{\int_{n}^{z_{r}}{{{dzf}_{RS}\left( {\overset{\_}{r}}_{b} \right)}^{1/2}{\sum\left( \overset{\_}{r} \right)^{1/2}}}}} = {f_{RS}\left( {\overset{\_}{r}}_{b} \right)}^{1/2}}}} & (6) \\ {\mspace{79mu}{{{z_{r} \approx {{f_{RS}\left( {\overset{\_}{r}}_{b} \right)}^{{- 1}/2}{\sum{\left( \overset{\_}{r} \right)^{{- 1}/2}{L_{p}.\mspace{79mu} x^{\text{?}}}}}}} = {x - {\overset{\_}{r}}_{b}}}\mspace{79mu}{{x^{\text{?}}/{\overset{\_}{r}}_{b}} ⪡ 1}\mspace{79mu}{{x^{\text{?}}/{\left. {\overset{\_}{r}}_{b} \right.\sim r_{S}}}/{\overset{\_}{r}}_{b}}\mspace{79mu}{r_{S}/{\overset{\_}{r}}_{b}}\mspace{79mu}{x^{\text{?}}/{\overset{\_}{r}}_{b}}}} & (7) \\ {\mspace{79mu}{g_{\nu,{\mu\text{?}}} = {{f_{RS}\left( {\overset{\_}{r}}_{b} \right)}\left( {1 + \frac{r_{S}x^{\text{?}}}{{\overset{\_}{r}}_{b}^{2}}} \right) \times}}} & (8) \\ {{{diag}\left( {{- 1},{\sum{\left( {\overset{\_}{r}}_{b} \right)\left( {1 - \frac{2r_{S}x^{\text{?}}}{{\overset{\_}{r}}_{b}^{2}}} \right)\text{?}\text{indicates text missing or illegible when filed}}}} \right.}} & (9) \end{matrix}$

The proper length in the vertical orientation is

$\begin{matrix} {\mspace{79mu}{L_{p} = {\int_{0}^{x\text{?}}{{dx}^{\text{?}}\left( {g_{\nu,11}(x)} \right)}^{1/2}}}} & (10) \\ {\mspace{110mu}{{= {{f_{RB}\left( {\overset{\_}{r}}_{b} \right)}^{1/2}{\sum{\left( {\overset{\_}{r}}_{b} \right)^{1/2}{x_{t}^{\text{?}}\left( {1 - \frac{r_{S}x_{t}^{\text{?}}}{4{\overset{\_}{r}}_{b}^{2}}} \right)}}}}},{\text{?}\text{indicates text missing or illegible when filed}}}} & (11) \end{matrix}$

which can be inverted as

$\begin{matrix} {\mspace{79mu}{x_{t}^{\text{?}} = {{f_{RS}\left( {\overset{\_}{r}}_{b} \right)}^{{- 1}/2}{\sum{\left( {\overset{\_}{r}}_{b} \right)^{{- 1}/2}{{L_{p}\left( {1 + \frac{r_{S}L_{p}}{4{\overset{\_}{r}}_{b}^{2}}} \right)}.\text{?}}\text{indicates text missing or illegible when filed}}}}}} & (12) \end{matrix}$

In the next Section, we will introduce the description of a BEC in the spacetime described by the metric (3) and we will give the conditions under which the phonon field in the BEC can be described completely analogously to the light field confined in an optical cavity placed in the same spacetime.

Section III. 1-Dimensional Relativistic BEC in the Gravitational Field of the Earth

A BEC can be effectively described by the Gross-Pitaevskii equation (see L. P. Pitaevskii, Soviet Physics JETP-USSR 13, 451 (1961) and E. P. Gross, Il Nuovo Cimento 20, 454 (1961)). This accounts for the collective behaviour of the bulk of the atoms that compose the BEC. However, small perturbations of the BEC, known as phonons, appear as quantised massless particles that propagate on an effective curved background metric. The kinematics of the phonons propagating through the bulk are described by a Klein-Gordon equation (21, below), where the metric g is replaced by the effective metric g with components

$\begin{matrix} {{\mathfrak{g}}_{\mu\nu} = {\frac{c}{c_{s}}{{\rho\left\lbrack {g_{\mu\nu} + {\left( {1 - \frac{c_{s}^{2}}{c^{2}}} \right)\frac{u_{\mu}u_{\nu}}{{u_{\alpha}u^{\alpha}}}}} \right\rbrack}.}}} & (13) \end{matrix}$

Here, the speed c_(s) is defined by

$\begin{matrix} {{c_{s}^{2}\text{:=}\mspace{14mu}\frac{c^{2}c_{0}^{2}}{{{u_{\alpha}u^{\alpha}}} + c_{0}^{2}}},{where}} & (14) \\ {c_{0}^{2}\text{:=}\mspace{14mu}\frac{h^{2}}{2m^{2}}{\lambda\rho}} & (15) \end{matrix}$

and the four-vector u_(μ), is the flow associated with the phase of the wave function of the BEC bulk (see C. Barceló, S. Liberati, and M. Visser, Classical and Quantum Gravity 18, 1137 (2001)), the speed of sound in the BEC is c₀ and the density ρ of the BEC may depend on space and time coordinates. For the BECs that we are considering here, c₀ ²/c² is of the order of 10⁻²⁰. Therefore, we can approximate c_(s)=c₀.

An important feature of BECs is that the background metric g can be experimentally tuned to replicate dynamical spacetimes. This can be done experimentally by tuning accordingly the density ρ, the flows u_(μ), and the speed of sound c₀. This versatile nature of BECs provides us with a tool for probing dynamics and properties of the “real” background spacetime when first attempting to study non-Minkowski background metrics.

The density, the flow and the speed of sound in the BEC cannot be chosen freely and are constrained. By considering two body collisions only, one can model a BEC in a curved spacetime by considering a charged, self interacting scalar field governed by the relativistic, non-linear field Eq. (see the generalisation of Eq. (10) in S. Fagnocchi, S. Finazzi, S. Liberati, M. Kormos, and A. Trombettoni, New J. Phys. 12, 095012 (2010) to curved spacetimes)

$\begin{matrix} {{{{\bullet\hat{\phi}} - {\left( {\frac{m^{2}c^{2}}{h^{2}} + V} \right)\hat{\phi}} - {\lambda{\hat{\phi}}^{2}\hat{\phi}}} = 0.}\left. \hat{\phi}\rightarrow\phi \right.{\hat{\phi} = {\sqrt{\rho}ɛ^{i\; 0}}}{u_{\mu} = {\frac{h}{m}{\partial_{\mu}\theta}}}} & (16) \\ {{0 = {\partial_{\mu}\left( \sqrt{{- {\det(g)}}g^{\mu\nu}\rho\; u_{\nu}} \right)}}{{u}^{2} = {{{- g^{\mu\nu}}u_{\mu}u_{\nu}} = {c^{2} + {{\frac{h^{2}}{m^{2}}\left\lbrack {V + {\lambda\rho} - \frac{\bullet\sqrt{\rho}}{\sqrt{\rho}}} \right\rbrack}.}}}}} & (17) \end{matrix}$

These equations connect the density ρ and the flows u_(μ).

As a first approach to the problem of analysing a trapped BEC placed in Schwarzschild spacetime, we assume that the external potential V allows for a constant density ρ and vanishing spatial components of u_(μ). As explained in Section II, we consider an effective 1+1-dimensional situation. To this end, we assume that the extension of the bulk is much larger in one direction than in the other two directions. The solution to Eqs. (17) in the elongated direction, when the cavity is horizontal and the trapping potential is uniform with infinite walls, i.e. V=0 inside the box and V=∞ outside of the box, becomes

u ₀ ² =f _(RS)( r _(b))(c ²+2c ₀ ²)≈c ² f _(RS)( r _(b)).  (18)

When considering the vertical orientation of the effective 1+1-dimensional BEC, the same expression for the flow component u₀ as in (18) can be obtained if a potential of the form V=V₀+ΔV is added, where V₀ is the uniform trapping potential and ΔV the linear potential given as

$\begin{matrix} {\mspace{79mu}{{\Delta\; V} = {{{- \frac{m^{2}}{\hslash^{2}}}\frac{r_{S}x^{\text{?}}}{{\overset{\_}{r}}_{b}^{2}}u_{0}^{2}} \approx {\frac{m^{2}c^{2}}{\hslash^{2}}{\frac{r_{S}x^{\text{?}}}{{\overset{\_}{r}}_{b}^{2}}.\text{?}}\text{indicates text missing or illegible when filed}}}}} & (19) \end{matrix}$

The potential ΔV cancels the effect of the gravitational field on the bulk up to second order in r_(S)/r _(b) and x/r _(b) and first order in their product. The cancellation of the gravitational field is common in experiments with BECs. It is done, for example, by applying a magnetic field of the right strength (see for instance A. L. Gaunt, T. F. Schmidutz, I. Gotlibovych, R. P. Smith, and Z. Hadzibabic, Phys. Rev. Lett. 110, 200406 (2013)). In order to obtain the potential ΔV, the experimenter has to tune the gradient of the linear potential, dΔV/dx, such that the density is homogeneous. However, dΔV/dx can also be fixed to a value that differs slightly from the value necessary to cancel the full effect of the gravitational field. This difference can lead to a density perturbation which has an effect on the phonon field that can, in principle, be measured. Let us assume that ΔV₀ is given by Eq. (19) but r_(S)/r _(b) ² replaced by some default value adjustments r_(S,0)/r _(b,0) ². Then, the corresponding relative density perturbation δ_(ρ) defined as ρ=ρ_(o)(1+δ_(ρ)) is given by

$\begin{matrix} {\mspace{76mu}{{{\delta_{\rho} = {{- \frac{1}{2}}\frac{c^{2}}{c_{\text{?}0}^{2}}\left( {\frac{r_{S}}{{\overset{\_}{r}}_{b}^{2}} - \frac{r_{S,0}}{{\overset{\_}{r}}_{b,0}^{2}}} \right)\left( {x^{\text{?}} - \frac{{x_{\text{?}}}_{t}^{\text{?}}}{2}} \right)}},\mspace{79mu}{{f_{0}^{x\text{?}}{dx}^{\text{?}}\delta_{p}} = 0}}{\text{?}\text{indicates text missing or illegible when filed}}}} & (20) \end{matrix}$

Section IV. Quantization of the Fields

In this section, we quantise the phonon/light field ψ in the BEC/cavity. The phonon/photon field obeys the massless Klein-Gordon equation (see N. D. Birrell and P. C. W. Davies, Quantum fields in curved space (Cambridge: University Press, 1982)

$\begin{matrix} {{{{\bullet\psi} = {\frac{1}{\sqrt{- g}}{\partial_{\mu}\left( {\sqrt{- g}g^{\mu\nu}\partial_{\nu}} \right)}}},{g\text{:=}\mspace{14mu}{\det\left( g_{\mu\nu} \right)}}}{g_{\mu\nu} = g_{\mu\nu}}{g_{\mu\nu} = {\mathfrak{g}}_{\mu\nu}}g_{\mu\nu}} & (21) \end{matrix}$

A. Horizontal Orientation

$\begin{matrix} {\mspace{70mu}{{{\left( {x^{0},\overset{\sim}{z}} \right)\mspace{14mu}{with}{\;\mspace{11mu}}\overset{\sim}{z}} = {\sum{\left( {\overset{\_}{r}}_{b} \right)^{{- 1}/2}z}}}\mspace{70mu}{{\overset{\sim}{g}}_{h,{\mu\nu}} = {{f_{RS}\left( {\overset{\_}{r}}_{b} \right)}{{{diag}\left( {{- 1},1} \right)}.\mspace{70mu}{\psi_{h}\left( {x^{0},{\overset{\sim}{z}}_{r}} \right)}}}}}} & (22) \\ {\mspace{70mu}{{{{\phi_{h,k}\left( {x^{0},\overset{\sim}{z}} \right)}\text{:=}\frac{e^{{- i}\;\omega_{h,k}t}}{\sqrt{\pi\; k}}{\sin\left( {\frac{w_{h,k}}{\text{?}}\overset{\sim}{z}} \right)}},\mspace{70mu}{\omega_{h,k}^{cav} = {\frac{c\;\pi\; k}{\text{?}} = {\frac{c\;\pi\; k}{L_{p}}{f_{RS}\left( {\overset{\_}{r}}_{b} \right)}^{1/2}}}}}{\text{?}\text{indicates text missing or illegible when filed}}}} & (23) \end{matrix}$

The solutions (23) are normalised with respect to the Klein-Gordon inner product

$\begin{matrix} {{\left( {\phi_{h},{\hat{\phi}}_{h}} \right)_{{KG},h}\text{:=}\mspace{14mu} i{\int_{0}^{2r}{d\overset{\sim}{z}{f\left( r_{b} \right)}^{1/2}\left( {{\phi_{h}^{*}{\partial_{0}{\overset{\sim}{\phi}}_{h}}} - {\phi_{h}^{*}{\partial_{0}\phi_{h}^{*}}}} \right)}}}{{{\psi_{h}\left( {x^{0},\overset{\sim}{z}} \right)} = {\sum\limits_{k}\left\{ {{{\hat{a}}_{h,k}{\phi_{h,k}\left( {x^{0},\overset{\sim}{z}} \right)}} + {{\hat{a}}_{h,k}^{\dagger}{\phi_{h,k}^{*}\left( {x^{0},\overset{\sim}{z}} \right)}}} \right\}}},}} & \; \\ {{\hat{a}}_{h,k},{{{\hat{a}}_{h,k}^{\dagger}\left\lbrack {{\hat{a}}_{h,k},{\hat{a}}_{h,k^{\prime}}^{\dagger}} \right\rbrack} = \delta_{k,h^{\prime}}}} & (24) \end{matrix}$

For the phonon field in the BEC, we have to consider the acoustic metric g from Eq. (13). Since the BEC is effectively 1-dimensional, we find, in the coordinates (c_(s)t z),

$\begin{matrix} {\mspace{70mu}{{\mathfrak{g}}_{h,{\rho\nu}} = {{\rho\frac{c}{c_{s}}{f_{RS}\left( {\overset{\_}{r}}_{b} \right)}{\left( {{- 1},{\sum\left( {\overset{\_}{r}}_{b} \right)}} \right).\mspace{70mu} w_{h,k}}} = {\frac{c_{s}\pi\; k}{L_{p}}{f_{RS}\left( {\overset{\_}{r}}_{b} \right)}^{1/2}}}}} & (25) \\ {\mspace{70mu}{{w_{h,k} = {\frac{\text{?}\;\pi\; k}{L_{p}}{f_{RS}\left( {\overset{\_}{r}}_{b} \right)}^{1/2}}},{\text{?}\text{indicates text missing or illegible when filed}}}} & (26) \end{matrix}$

where ν=c for the light field and ν=c_(s) for the phonon field.

B. Vertical Orientation

$\begin{matrix} {\mspace{76mu}{\left( {x^{0},\overset{\sim}{x}} \right)\mspace{70mu}{{\overset{\sim}{x} = {{\sum{\left( {\overset{\_}{r}}_{b} \right)^{1/2}{\int_{0}^{x^{\text{?}}}{{dx}^{''}\left( {1 - \frac{2r_{S}x^{\text{?}}}{{\overset{\_}{r}}_{b}^{2}}} \right)}^{1/2}}}}\mspace{79mu} \approx {\sum{\left( {\overset{\_}{r}}_{b} \right)^{1/2}{x^{f}\left( {1 - \frac{r_{S}x^{\prime}}{2{\overset{\_}{r}}_{b}^{2}}} \right)}}}}},\mspace{70mu}{x^{\text{?}} = {\sum{\left( {\overset{\_}{r}}_{b} \right)^{{- 1}/2}{\hat{x}\left( {1 + {r_{S}{x^{\text{?}}/2}{\overset{\sim}{r}}_{b}^{2}}} \right)}}}},}}} & (27) \\ {\mspace{70mu}{{{\hat{g}}_{\nu,{\mu\upsilon}} = {{f_{RS}\left( {\overset{\sim}{r}}_{b} \right)}\left( {1 + \frac{r_{S}\hat{x}}{{\overset{\_}{r}}_{b}^{2}}} \right){{diag}\left( {{- 1},1} \right)}}}\mspace{70mu}{r_{S}/{\overset{\_}{r}}_{b}^{2}}\mspace{70mu}{x^{\text{?}}/{\overset{\_}{r}}_{b}^{2}}}} & (28) \\ {\mspace{70mu}{{{\psi_{v}\left( {x^{0},{\overset{\sim}{x}}_{b}} \right)} = {{\psi_{v}\left( {x^{0},{\overset{\sim}{x}}_{t}} \right)} = 0}}\mspace{70mu}{{{{\hat{\psi}}_{v}\left( {x^{0},\overset{\_}{x}} \right)} = {\sum\limits_{k}\;{\left\{ {{{\hat{a}}_{v,k}{\phi_{v,k}\left( {x^{0},\overset{\_}{x}} \right)}} + {{\hat{a}}_{h,k}^{\dagger}{\phi_{v,k}^{*}\left( {x^{0},\overset{\_}{x}} \right)}}} \right\}.\mspace{70mu}{\hat{a}}_{v,k}}}},{{{\hat{a}}_{v,k}^{\dagger}\mspace{70mu}\left\lbrack {{\hat{a}}_{v,k},{\hat{a}}_{v,k^{\prime}}^{\dagger}} \right\rbrack} = \delta_{k,k^{\prime}}}}}} & (29) \\ {\mspace{70mu}{{{{\phi_{v,k}\left( {x^{0},\overset{\sim}{x}} \right)}\text{:=}\mspace{14mu}\frac{e^{{- i}\;\omega_{v,k}t}}{\sqrt{\pi\; k}}{\sin\left( {\frac{\omega_{v,k}}{c}\overset{\_}{x}} \right)}},\mspace{70mu}{\omega_{v,k}^{cav}\text{:=}\mspace{14mu}\frac{c\;\pi\; k}{{\overset{\sim}{x}}_{t}}}}\mspace{70mu}{\overset{\sim}{x}}_{t}}} & (30) \\ {\mspace{70mu}{{\overset{\sim}{x}}_{t} \approx {{f_{RS}\left( {\overset{\_}{r}}_{b} \right)}^{{- 1}/2}{{L_{p}\left( {1 - \frac{r_{S}L_{p}}{4{\overset{\sim}{r}}_{b}^{2}}} \right)}.\text{?}}\text{indicates text missing or illegible when filed}}}} & (31) \end{matrix}$

The phonon field can be treated similarly to the horizontal case if the external linear potential is exactly canceling the effect of the gravitational field, i.e. if δ_(ρ)=0. However, to fix the potential, a precise measurement of the gravitational field would be necessary in the first place which would reduce the utility of our setup for measurement significantly. Instead, the external potential can be fixed to some value close to the cancellation and the effect of the difference can be used for the measurement of the gravitational field itself. As explained in Section III, there are no spatial velocity flows, and we find for the acoustic metric due to a time-independent density perturbation

$\begin{matrix} {\mspace{76mu}{{\mathfrak{g}}_{v,{\mu\upsilon}} \approx {\rho_{0}\frac{c}{c_{s\; 0}}\left( {1 + \delta_{p}} \right)^{1/2}{\begin{pmatrix} {\frac{c_{s\; 0}^{2}}{c^{2}}\left( {1 + \delta_{\rho}} \right)^{1/2}g_{v,00}} & 0 \\ 0 & g_{v,11} \end{pmatrix}.}}}} & (32) \\ {\mspace{76mu}{{x^{0} = {c_{s\; 0}t}}\mspace{76mu}{{\overset{\_}{\mathfrak{g}}\;}_{v,{\mu\upsilon}} = {\rho_{0}\frac{c}{c_{s\; 0}}\left( {1 + \delta_{p}} \right)^{1/2}\begin{pmatrix} {\left( {1 + \delta_{p}} \right)^{1/2}g_{v,00}} & 0 \\ 0 & g_{v,11} \end{pmatrix}}}\mspace{76mu}\left( {x^{0},x^{\text{?}}} \right)\mspace{76mu}{x^{\text{?}} = {x = {\overset{\sim}{r}}_{b}}}\mspace{76mu}{{x^{\text{?}}/{\overset{\_}{r}}_{b}} ⪡ 1}\mspace{76mu}{{x^{\text{?}}/{\left. {\overset{\_}{r}}_{b} \right.\sim r_{S}}}/{\overset{\sim}{r}}_{b}}\mspace{76mu}{r_{S}/{\overset{\_}{r}}_{b}}\mspace{76mu}{x^{\text{?}}/{\overset{\_}{r}}_{b}}}} & (33) \\ {\mspace{76mu}{{\overset{\_}{\mathfrak{g}}}_{v,00} = {- {P\left( x^{\text{?}} \right)}}}} & (34) \\ {\mspace{76mu}{{{\overset{\_}{\mathfrak{g}}}_{v,11} = {{- {P\left( x^{\text{?}} \right)}}\left( {1 - \delta_{p} + \delta_{p}^{2}} \right){\sum{\left( {\overset{\_}{r}}_{b} \right)\left( {1 - \frac{2r_{S}x^{\text{?}}}{{\overset{\_}{r}}_{b}^{2}}} \right)}}}}\mspace{76mu}{where}\mspace{76mu}{{P\left( x^{\text{?}} \right)} = {\rho_{0}\frac{c}{c_{s\; 0}}\left( {1 + \frac{3\delta_{p}}{2} + \frac{3\delta_{p}}{8}} \right){f_{RS}\left( {\overset{\_}{r}}_{b} \right)}{\left( {1 + \frac{r_{S}x^{\text{?}}}{{\overset{\_}{r}}_{b}^{2}}} \right).}}}}} & (35) \\ {\mspace{76mu}\left( {x^{0},\overset{\sim}{x}} \right)} & (36) \\ {\mspace{76mu}{\overset{\sim}{x}\text{:=}\mspace{14mu}{\int_{0}^{x^{\text{?}}}{{dx}^{''}{\sqrt{\frac{{\overset{\_}{\mathfrak{g}}}_{v,11}\left( x^{''} \right)}{{{\overset{\_}{\mathfrak{g}}}_{v,00}\left( x^{''} \right)}}}.\text{?}}\text{indicates text missing or illegible when filed}}}}} & (37) \end{matrix}$

In these coordinates, the acoustic metric becomes

$\begin{matrix} {\mspace{76mu}{{{\overset{\sim}{\mathfrak{g}}}_{v,{\mu\upsilon}} = {{\hat{P}\left( \overset{\sim}{x} \right)}{{diag}\left( {{- 1},1} \right)}}},\mspace{76mu}{{\overset{\sim}{P}\left( \overset{\sim}{x} \right)} = {P\left( {x^{\text{?}}\left( \overset{\_}{x} \right)} \right)}}}} & (38) \\ {\mspace{76mu}{{{{\phi_{v,k}\left( {x^{0},\overset{\_}{x}} \right)} = {\frac{e - {i\;\omega_{\text{?},k}t}}{\sqrt{\pi\; k}}{\sin\left( \frac{\omega_{v,k}}{c_{s\; 0}} \right)}}},\mspace{76mu}{\omega_{v,k}^{BEC} = \frac{c_{s\; 0}\pi\; k}{{\hat{x}}_{t}}}}\mspace{76mu}{\hat{x}}_{t}\mspace{76mu}{\hat{x}}_{t}\mspace{76mu}{{\hat{x}}_{t} = {{\sum{\left( {\overset{\_}{r}}_{b} \right)^{1/2}{\int_{0}^{x\text{?}}{{{dx}^{\text{?}}\left( {1 - \frac{\delta_{\rho}}{2} + \frac{3\delta_{\rho}}{8}} \right)}\left( {1 - \frac{r_{S}x^{\text{?}}}{{\overset{\_}{r}}_{b}^{2}}} \right)}}}}\mspace{85mu} = {\sum{\left( {\overset{\_}{r}}_{b} \right)^{1/2}{\int_{0}^{x\text{?}}{{dx}^{\text{?}}\left( {1 - \frac{r_{S}x^{\text{?}}}{{\overset{\_}{r}}_{b}^{2}} + {\frac{r_{S}x^{\text{?}}}{{\overset{\_}{r}}_{b}^{2}}\frac{\delta_{\rho}}{2}} + \mspace{115mu}{\left( {1 - \frac{r_{S}x^{\text{?}}}{{\overset{\sim}{r}}_{b}^{2}}} \right)\frac{3\delta_{\rho}}{8}}} \right)}}}}}}}} & (39) \\ {\mspace{76mu}{\delta_{\rho}\mspace{76mu}{r_{S}{L_{p}/{\overset{\_}{r}}_{b}^{2}}}\mspace{76mu}\delta_{\rho}\mspace{76mu}{r_{S}{L_{p}/{\overset{\_}{r}}_{b}^{2}}}\mspace{76mu}\delta_{\rho}\mspace{76mu}{r_{S}{L_{p}/{\overset{\_}{r}}_{b}^{2}}}\mspace{76mu}{c^{2}/c_{s\; 0}^{2}}\mspace{76mu}\delta_{\rho}\mspace{76mu}{r_{S}{L_{p}/r_{b}^{2}}}\mspace{76mu}{r_{S}/{\overset{\_}{r}}_{b}^{2}}{\text{?}\text{indicates text missing or illegible when filed}}}} & (40) \end{matrix}$

r_(S)/r _(b) ² and δ_(ρ) ² in Eq. (40) in the following. From Eq. (20), we see that δ_(p)=δ′_(ρ)(x′−x′_(t)/2), where δ′_(ρ)=dδ_(ρ)/dx′. We find

$\begin{matrix} \begin{matrix} {{\overset{\sim}{x}}_{t} \approx {\sum{\left( {\overset{\_}{r}}_{b} \right)^{1/2}{\int_{0}^{x_{t}^{\prime}}{{dx}^{\prime}\left( {1 - \frac{r_{S}x^{\prime}}{{\overset{\_}{r}}_{b}^{2}} + {\frac{r_{S}}{{\overset{\_}{r}}_{b}^{2}}\frac{\delta_{\rho}^{\prime}}{2}\left( {{x^{\prime}}^{2} - \frac{x^{\prime}x_{t}^{\prime}}{2}} \right)} +} \right.}}}}} \\ \left. {\frac{3\delta_{\rho}^{\prime}}{8}\left( {{x^{\prime}}^{2} - {x^{\prime}x_{t}^{\prime}} + \frac{{x^{\prime}}_{t}^{2}}{4}} \right)} \right) \\ {= {\sum{\left( {\overset{\_}{r}}_{b} \right)^{1/2}{x_{t}^{\prime}\left( {1 - \frac{r_{S}x_{t}^{\prime}}{2{\overset{\_}{r}}_{b}^{2}} + {\frac{r_{S}}{{\overset{\_}{r}}_{b}^{2}}\frac{\delta_{\rho}^{\prime}}{24}{x^{\prime}}_{t}^{2}} + {\frac{\delta_{\rho}^{\prime 2}}{64}{x^{\prime}}_{t}^{2}}} \right)}}}} \end{matrix} & (41) \end{matrix}$

With Eq. (12), we find

$\begin{matrix} {{{\overset{\sim}{x}}_{t} \approx {{f_{RS}\left( {\overset{\_}{r}}_{b} \right)}^{{- 1}/2}{{{L_{p}\left( {1 - \frac{r_{S}L_{p}}{4{\overset{\_}{r}}_{b}^{2}} + {\left( {{\frac{r_{S}}{{\overset{\_}{r}}_{b}^{2}}\frac{\delta_{\rho}^{\prime}}{24}} + \frac{\delta_{\rho}^{\prime}}{64}} \right)L_{p}^{2}}} \right)}.c^{2}}/c_{s\; 0}^{2}}}}{{r_{S}/{\overset{\_}{r}}_{b}^{2}}\delta_{\rho}^{\prime}}{{\overset{\sim}{x}}_{t} \approx {{f_{RS}\left( {\overset{\_}{r}}_{b} \right)}^{{- 1}/2}{L_{p}\left( {1 - \frac{r_{S}L_{p}}{4{\overset{\_}{r}}_{b}^{2}} + \frac{{\delta^{\prime}}_{\rho}^{2}L_{p}^{2}}{64}} \right)}}}{{\delta_{V} \sim 10^{- 6}},{{c^{2}/c_{s\; 0}^{2}} \sim {10^{20}L_{p}} \sim 10^{- 4}}}{{{\delta^{\prime}}_{p}^{2}{L_{p}^{2}/64}} \sim 10^{- 6}}{{r_{S}{L_{p}/4}{\overset{\_}{r}}_{b}^{2}} \sim 10^{- 21}}\delta_{p}{r_{S}{L_{p}/{\overset{\_}{r}}_{b}^{2}}}} & (42) \end{matrix}$

C. Frequency Shifts

$\begin{matrix} {{{r\overset{\_}{r}}r = {\left( {1 + {{r_{S}/4}\overset{\_}{r}}} \right)^{2}\overset{\_}{r}}}{{f_{RS}\left( \overset{\_}{r} \right)} = {{1 - {r_{S}/r}} = {{f(r)}\;{r_{S}/r}}}}{\omega_{h,k}^{cav} = {\frac{c\;\pi\; k}{L_{p}}{f\left( r_{b} \right)}^{1/2}\mspace{14mu}{and}}}} & (43) \\ {\omega_{v,k}^{cav} = {\omega_{h,k}^{cav}\left( {1 + \frac{r_{S}L_{p}}{4\; r_{b}^{2}}} \right)}} & (44) \end{matrix}$

for the light cavity and

$\begin{matrix} {\omega_{h,k}^{BEC} = {\frac{{c\;}_{s\; 0}\pi\; k}{L_{p}}{f\left( r_{b} \right)}^{1/2}\mspace{14mu}{and}}} & (45) \\ {\omega_{v,k}^{BEC} = {\omega_{h,k}^{BEC}\left( {1 + \frac{r_{S}L_{p}}{4\; r_{b}^{2}} - \frac{{\delta^{\prime}}_{\rho}^{2}L_{p}^{2}}{64}} \right)}} & (46) \end{matrix}$

for the BEC and where c_(s)=c₈₀ for the horizontal cavity. The factor f(r_(b))¹/² represents the usual gravitational red shift that is experienced when photons are sent between two points at different height in the gravitational field of the earth (see R. V. Pound and G. A. Rebka, Phys. Rev. Lett. 3, 439 (1959)). Note that, due to the red shift, both frequencies become zero at the Schwarzschild radius, which occurs when approaching a black hole. Hence, this limit has to be taken with care. For our considerations about the gravitational field of the Earth, this is not of importance since every point outside of the surface of the Earth is far from its Schwarzschild radius.

In the limit L_(p)→0, we have ω_(v,k)→ω_(h,k). Therefore on the one hand, we conclude that the root of the frequency is the finite spatial extension of the BEC/cavity. In particular, the frequencies ω_(v,k) for a BEC/cavity with its lower end at r_(b) are equivalent to the frequencies of the field modes in the horizontal BEC/cavity at r_(e)=r_(b)+L_(P)/2 in second order in L_(P)/r_(b) and in first order in their product. Hence, there would be no difference between the two systems in this order if the horizontal cavity would be located in the middle between the ends of the vertical cavity. On the other hand, the frequency difference between ω_(v,k) and ω_(h,k) for both cavity and BEC, is due to the proper acceleration of the whole system that is necessary to keep it at a fixed Schwarzschild coordinate r. In particular, ω_(v,k)=ω_(h,k) only at r_(b)=∞ that is, asymptotically far from the earth.

For the derivation Eq. (44), we assumed that the proper length of the BEC/cavity is kept fixed by taking a perfectly rigid rod as reference. However, all realistic materials are deformable. In the next section, we will derive the effect of the deformation of the reference rod on the frequency spectrum of the fields in the gravitational field of the earth.

Section V. Deformation of the Cavity

We assume that the cavity/trap is supported by a rod of constant but negligible width. Furthermore, we assume that the material properties are constant along the rod. To keep the cavity/trap at a constant distance from the centre of the Earth, the whole system has to be accelerated upwards. We assume that this is achieved by supporting the rod from below at r_(b) (not by suspending it from above).

$\begin{matrix} {\mspace{79mu}{{{\left( {\text{?},\overset{\_}{x}} \right)\mspace{79mu}{{r_{S}/{\overset{\_}{r}}_{b}}{\overset{\_}{x}/{\overset{\_}{r}}_{b}}}\mspace{79mu}\frac{d^{2}}{d\;\tau^{2}}{\xi^{\mu}(\tau)}} = {{- {\Gamma_{\alpha\;\beta}^{\mu}\left( {\xi(\tau)} \right)}}{{\overset{.}{\xi}}^{\alpha}(\tau)}{{\overset{.}{\xi}}^{\beta}(\tau)}}},{\text{?}\text{indicates text missing or illegible when filed}}}} & (47) \end{matrix}$

where the dot denotes the derivative with respect to the proper time T of the segment and the connection coefficients are given as

$\begin{matrix} {\Gamma_{v\;\rho}^{\mu} = {{\frac{1}{2}{{{\overset{\sim}{g}}_{v}^{\mu\;\sigma}\left( {{\partial_{v}{\overset{\sim}{g}}_{v,{\sigma\;\rho}}} + {\partial_{\rho}{\overset{\sim}{g}}_{v,{\sigma\; v}}} - {\partial_{\rho}{\overset{\sim}{g}}_{v,{v\;\rho}}}} \right)}.{{\overset{.}{\xi}}^{\beta}(r)}}} = \left( {{()}^{{- \lambda}/2},0} \right)}} & (48) \\ {{a^{\overset{\sim}{x}} = {\frac{d^{2}\xi^{x}}{d\; t^{2}} = {{\left( {\overset{.}{\xi}}^{0} \right)^{- 2}\frac{d^{2}\xi^{x}}{d\;\tau^{2}}} = {- \frac{c^{2}r_{S}}{2{\overset{\_}{r}}_{b}^{2}}}}}}{\overset{\sim}{x} \in \left\lbrack {0,{\overset{\sim}{x}}_{t}} \right\rbrack}d{\overset{\sim}{x}\left( {\overset{\sim}{g}}_{v,11} \right)}^{1/2}d{\overset{\sim}{x}\left( {\overset{\sim}{g}}_{v,11} \right)}^{1/2}\rho\; d\overset{\sim}{x}} & (49) \\ {\begin{matrix} {{\sigma_{\overset{\sim}{x}}\left( \overset{\sim}{x} \right)} = {a^{{\overset{\sim}{x}}_{t}}{\int_{\overset{\sim}{x}}^{{\overset{\sim}{x}}_{t}}{d{\overset{\sim}{x}\left( {\overset{\sim}{g}}_{v,11} \right)}^{1/2}\rho}}}} \\ {= {{f_{RS}^{1/2}\left( {\overset{\_}{r}}_{b} \right)}\rho\;{{a^{\overset{\sim}{x}}\left( {{\overset{\sim}{x}}_{t} - \overset{\sim}{x}} \right)}.}}} \end{matrix}{{\left( \overset{\sim}{x} \right)} = {{\sigma_{\overset{\sim}{x}}\left( \overset{\sim}{x} \right)}/E}}r = {\left( {1 + {{r_{S}/4}\overset{\_}{r}}} \right)^{2}\overset{\_}{r}}} & (50) \\ {{\Delta\; L_{p}} = {\int_{0}^{{\overset{\sim}{x}}_{t}}{d{\overset{\sim}{x}\left( {\overset{\sim}{g}}_{v,11} \right)}^{1/2}{\left( \overset{\sim}{x} \right)}}}} & (51) \\ {= {{f_{RS}\left( {\overset{\_}{r}}_{b} \right)}\frac{\rho\; a^{\overset{\sim}{x}}{\overset{\sim}{x}}_{t}^{2}}{2\; E}}} & (52) \\ {{\approx {{- \frac{c^{2}\rho}{E}}\frac{r_{S}L_{p}^{2}}{4{\overset{\_}{r}}_{b}^{2}}} \approx {{- \frac{c^{2}\rho}{E}}\frac{r_{S}L_{p}^{2}}{4r_{b}^{2}}}}{E/\rho_{rod}}} & (53) \\ {{\Delta\; L_{p}} = {{- \frac{c^{2}}{c_{rod}^{2}}}{\frac{r_{S}L_{p}^{2}}{4r_{b}^{2}}.}}} & (54) \end{matrix}$

Another meaningful question arising is, will tidal effects on the horizontal rod take part in the sensitivity? Due to the sphericity of the planet, the horizontal rod will be an effective cantilevered rigid beam supported on z_(l)=0 and subjected to gravitational acceleration on every other point of the system, as depicted in FIG. 3. The classical displacement δx at each point z∈┌0, z_(r)┐ of the cavity is effectively given by equation (2.3) of V. Srivastava, H. Jones, and G. Greenwood, Proceedings of the Royal Society of London A: Mathematical, Physical and Engineering Sciences 462, 2863 (2006). By considering a constant gravitational force (up to second order in z_(r)/r_(b), z/r_(b) and first order in their product), the displacement is

$\begin{matrix} {{{\delta\; x} = {3\frac{c^{2}}{c_{rod}^{2}}\frac{r_{S}}{r_{b}^{2}}\frac{1}{h^{2}}\left( {\frac{z_{r}^{2}z^{2}}{2} - \frac{z_{r}z^{3}}{3} + \frac{z^{4}}{12}} \right)}},{L_{p} \sim z_{r}}} & (55) \\ \begin{matrix} {{d\; L_{p}^{tidal}} = \sqrt{{d\; z^{2}} + {d\left( {\delta\; x} \right)}^{2}}} \\ {\approx {\left( {1 + \left( \frac{d\left( {\delta\; x} \right)}{d\; z} \right)^{2}} \right){{dz}.}}} \end{matrix} & (56) \end{matrix}$

Although this is clearly not the length that the rod will have in reality, it will serve with as a heuristic argument for our purposes.

$\begin{matrix} {{\omega_{h,k}^{tidal}:={v\;\pi\;{k/L_{p}^{tidal}}}}{\omega_{h,k}^{tidal} \approx {{{\omega_{h,k}\left( {1 - {\frac{9}{28}\frac{c^{4}}{c_{rod}^{4}}\frac{r_{S}^{2}}{r_{b}^{4}}\frac{L_{p}^{6}}{h^{4}}}} \right)}.3}\; L_{p}{r_{S}/2}\; r_{b}^{2}}}{3L_{p}c^{2}{r_{S}/2}c_{rod}^{2}r_{b}^{2}}{{\delta^{\prime}}_{\rho}^{2}{L_{p}^{2}/64}}} & (57) \end{matrix}$

For our purposes we will assume that the rod is sufficiently thick for this effect to be negligible. In the case of a light cavity, the change of proper length L_(p) only affects the radar length, and we can write the resonant frequencies as

$\begin{matrix} {\omega_{v,k}^{cav} = {{\omega_{h,k}^{cav}\left( {1 + \Lambda} \right)}.{where}}} & (58) \\ {\Lambda:={\frac{1}{4}\frac{r_{S}}{r_{b}^{2}}\left( {\frac{c^{2}}{c_{rod}^{2}} + 1} \right){L_{p}.}}} & (59) \end{matrix}$

For an aluminium rod, c²/c² _(rod) is of the order 10⁸ and for a rod of carbyne (linear acetylenic carbon), the material with the highest known specific modulus, it would be of the order 10⁸. Therefore, in all realistic situations, c²/c² _(rod)+1 is equivalent to c²/c² _(rod) in equation (59).

In the case of a uniformly trapped BEC, the length change also affects the speed of sound c, since it is proportional to the square root of the density of the BEC which is proportional to the square root of the inverse of the proper length of the trap. We obtain that

$\begin{matrix} {\omega_{v,k}^{BEC} = {{\omega_{h,k}\left( {1 + {\frac{3}{2}\Lambda}} \right)}.}} & (60) \end{matrix}$

We write the general equation

$\begin{matrix} {\omega_{h,k}^{def} = {\frac{v\;\pi\; k}{L_{p}}{f\left( r_{b} \right)}^{1/2}\mspace{14mu}{and}}} & (61) \\ {\omega_{v,k}^{def} = {{\omega_{h,k}\left( {1 + \Theta} \right)}.}} & (62) \end{matrix}$

where v and Θ are defined for the light cavity and the BEC respectively as

$\begin{matrix} {{Light}\mspace{14mu}\left\{ {\begin{matrix} {v = c} \\ {\Theta = \Lambda} \end{matrix},{{BEC}\;\left\{ {\begin{matrix} {v = c_{s\; 0}} \\ {\Theta = {{\frac{3}{2}\Lambda} - \frac{{\delta^{\prime}}_{\rho}^{2}L_{p}^{2}}{64}}} \end{matrix}.} \right.}} \right.} & (63) \end{matrix}$

The time evolution of a single mode of the field gives rise to a phase which can be employed to measure time. In this sense, these modes correspond to extended quantum clocks. We have already noted that ω_(h,k)→ω_(v,k) for r_(b)→∞. Therefore, the extended clocks work like local clocks for the observer at an infinite radial distance. This occurs because, at an extremely large (infinite) distance from the surface of the Earth no proper acceleration is needed to keep the system at a fixed coordinate r. We will discuss further properties of our extended clocks in Section A (below). With two clocks at different height, the red shift and the Schwarzschild radius of the earth can be inferred. This will be part of Section A as well. However, also the frequency difference of modes of two differently oriented BECs/cavities due to proper acceleration can employed for a measurement of that acceleration. When the radius to the center of the earth is known, the Schwarzschild radius of the earth can be inferred. We will discuss this application in Sections VI and VII.

A. Measurements of Proper Time

Special relativity predicts that time is not an absolute concept. However, a meaningful notion of time can be constructed as the time measured by a point-like clock moving along the world-line of a physical (pointlike) observer. This time is known as proper time r and is defined as τ=−∫ds, where ds is the line element. In the gravitational field of the Earth, we have τ_(r)=f(r)^(1/2)t for an observer at a fixed radial position r. In the following, we want to calculate the principle error bound on the estimation of the proper time τ_(rb) by a measurement on a single mode of the quantum field in the extended clock built from an optical cavity or the phonon field in a BEC. To this end, we employ quantum metrology and the quantum Cramer-Rao bound which gives a lower bound for the absolute error of the estimation of a parameter λ from a measurement on a quantum system Δλ² given a specific initial state. We have for the lower bound

$\begin{matrix} \left. {\left. \mspace{79mu}{{{{{\Delta\;\lambda^{2}} \geq \frac{1}{{MF}\left( {\lambda,Q} \right)} \geq \frac{1}{{MH}(\lambda)}},\mspace{79mu}{{{{{\Delta\;\lambda}}/{\lambda }} \geq {{\delta(\lambda)}{\delta(\lambda)}}}:=\left( {{\lambda }\sqrt{{MH}(\lambda)}} \right)^{- 1}}}\mspace{79mu}{{\hat{D}}_{k}(\alpha)}{{\hat{S}}_{k}(r)}}❘0_{k}} \right\rangle\mspace{79mu}{{{\hat{S}}_{k}(r)} = {{{\exp\left\lbrack {\frac{r}{\text{?}}\left( {{\hat{a}}_{k}^{\dagger 2} - {\hat{a}}_{k}^{2}} \right)} \right\rbrack}r} \in {\mathbb{R}}}}\mspace{79mu}{{{\hat{D}}_{k}(\alpha)} = {{\exp\left\lbrack {{\alpha{\hat{a}}_{k}^{\dagger 2}} - {\alpha^{*}{\hat{a}}_{k}}} \right\rbrack}\mspace{85mu} ❘0_{k}}}} \right\rangle & (64) \\ {\mspace{79mu}{{{\delta_{k}^{SMSC}\left( \tau_{r_{b}} \right)} = \frac{1}{\sqrt{2\; M}{{\sqrt{{\sinh^{2}\left( {2r} \right)} + {2{\alpha }^{2}e^{2r}}}\tau_{r_{b}}\frac{d}{d\;\tau_{r_{b}}}\psi_{k}}}}},\mspace{79mu}{\psi_{h,k} = {{\omega_{h,k}t} = {\omega_{0,k}\tau_{r_{b}}}}}}} & (65) \\ {\mspace{79mu}{{\psi_{v,k} = {{\omega_{v,k}t} = {{\omega_{0,k}\left( {1 + \Theta} \right)}\tau_{r_{b}}}}}\mspace{85mu}{{\omega_{0,k}:=\frac{v\;\pi\; k}{L_{p}}},{\text{?}\text{indicates text missing or illegible when filed}}}}} & (66) \end{matrix}$

which is the k-th frequency of the field in the cavity when it is placed asymptotically far from the Earth. Employing Eq. (61) and Eq. (62), we obtain the relative error bound for proper time measurement with the SMSC state to first order in Λ, which reads

$\begin{matrix} {{\delta_{h,k}^{SMSC}\left( \tau_{r_{b}} \right)} = \frac{1}{\sqrt{2\; M}\sqrt{{\sinh^{2}\left( {2r} \right)} + {2{\alpha }^{2}e^{2r}{{\omega_{0,k}\tau_{r_{b}}}}}}}} & (67) \\ {{\delta_{v,k}^{SMSC}\left( \tau_{r_{b}} \right)} \approx {{\delta_{h,k}^{SMSC}\left( \tau_{r_{b}} \right)}{\left( {1 - \Theta} \right).}}} & (68) \end{matrix}$

Given our results (67) and (68), we can answer the following question: How well can the proper time τ_(rb) be measured using an extended clock (cavity) with its lower end at r_(b) in the vertical setup and the whole clock at r_(b) in the horizontal setup?

a _(g) =c ² r _(a)/(2r _(b) ²)

B. Estimation of Red Shift

$\begin{matrix} {\mspace{79mu}{{{r_{b}\tau_{b^{\prime}}}\mspace{79mu}{{r_{b^{\prime}}r_{b}}\mspace{79mu} = {{\left( {\text{?} - \text{?}} \right)/\text{?}}\;\lambda_{r_{b}}\lambda_{r_{b^{\prime}}}}}}\mspace{79mu}{r_{b}r_{b^{\prime}}}\mspace{79mu}{\lambda_{r_{b}} = {{2\;\pi\;{c/\omega_{r_{b}}}\text{?}} = {\omega_{r_{b}}\text{?}}}}\mspace{79mu}{z = {\frac{T_{r_{b^{\prime}}} - T_{r_{b}}}{T_{r_{b}}} = {\frac{T_{r_{b^{\prime}}}}{T_{r_{b}}} - 1.}}}{\text{?}\text{indicates text missing or illegible when filed}}}} & (69) \end{matrix}$

From standard uncertainty propagation, we know that the relative error bound behaves as

$\begin{matrix} {{{\delta(z)} = {\sqrt{{\left( {\frac{\tau\text{?}}{z}\frac{\partial z}{{\partial\tau}\text{?}}} \right)^{2}\delta\text{?}\left( {\tau\text{?}} \right)^{2}} + {\left( {\frac{\tau_{r_{b}}}{z}\frac{\partial z}{{\partial\tau}\text{?}}} \right)^{2}\delta\text{?}\left( \tau_{r_{b}} \right)^{2}}} = {{\frac{{f\left( \text{?} \right)}^{1/2}}{{{f\left( \text{?} \right)}^{1/2} - {f\left( \text{?} \right)}^{1/2}}}{\sqrt{{\delta\text{?}\left( {\tau\text{?}} \right)^{2}} + {\delta\text{?}\left( {\tau_{r_{b}}\text{?}} \right)^{2}}} \cdot \text{?}}} = \text{?}}}}\mspace{79mu}{\text{?} = r_{b}}\mspace{79mu}\left. {r\text{?}}\rightarrow\infty \right.\mspace{79mu}\left. \text{?}\rightarrow\infty \right.\mspace{85mu}\left. \left( {{f\left( \text{?} \right)}\text{/}{f\left( {r\text{?}} \right)}} \right)^{1\text{/}2}\rightarrow{f\left( r_{b} \right)}^{1\text{/}2} \right.\mspace{79mu}\tau_{b}\mspace{79mu}{\omega_{0,k} = {\text{?}/L_{p}}}\mspace{20mu}{{\delta\left( z_{\infty} \right)} = {\frac{1}{{1 - {f\left( r_{b} \right)}^{1\text{/}2}}}{\sqrt{{\delta_{\infty}\left( \tau_{\infty} \right)}^{2} + {\delta_{r_{b}}\left( \tau_{r_{b}} \right)}^{2}}.\text{?}}\text{indicates text missing or illegible when filed}}}} & (70) \end{matrix}$

In the next Section, we discuss the estimation of the Schwarzschild radius of the Earth using the extended clocks.

C. Estimation of the Schwarzschild Radius r_(s)

$\begin{matrix} {\mspace{79mu}{{{{f\left( r_{b} \right)} = \left( {1 - {r_{S}/r_{b}}} \right)}\mspace{20mu} r_{S}\mspace{20mu}\tau_{b^{\prime}}\mspace{20mu}\tau_{b}\mspace{20mu}\tau_{b^{\prime}}\mspace{20mu}\text{?}}\mspace{20mu}\tau_{r_{b} - \tau_{r_{b^{\prime}}}}\mspace{20mu}\tau_{b}\mspace{20mu}{\tau\text{?}}\mspace{20mu}{z_{m}\tau_{{{r_{V}/r}\text{?}} - 1}}{{\delta\left( r_{S} \right)} = {{\frac{z}{r_{S}}{\frac{\partial r_{S}}{\partial z}}{\delta(z)}} = {{\frac{z}{r_{S}}{\frac{\partial z}{\partial r_{S}}}^{- 1}{\delta(z)}} = {{\frac{2r_{b}}{r_{S}}{f\left( r_{b} \right)}{\sqrt{\frac{{\delta_{r_{b^{\prime}}}\left( \tau_{r_{b^{\prime}}} \right)}^{2} + {\delta_{r_{b}}\left( \tau_{r_{b}} \right)}^{2}}{{1 - {\frac{r_{b}}{r_{b^{\prime}}}\frac{f\left( r_{b} \right)}{f\left( r_{b^{\prime}} \right)}}}}} \cdot \text{?}}} = {{\tau_{b}\mspace{11mu} r_{b^{\prime}}} = r_{b}}}}}}{\text{?}\text{indicates text missing or illegible when filed}}}} & (71) \end{matrix}$

the standard general relativity tenet that locally, space-time is flat. Let us consider the particular case of the emitter at an infinite radial distance from the surface of the Earth r_(b′)→∞. From Eq. (71), we obtain

$\begin{matrix} {{\delta_{k}^{{SMSV},\infty}\left( r_{S} \right)} = {2\left( {\frac{r_{b}}{r_{S}} - 1} \right)\sqrt{{\delta\left( \tau_{r_{\infty}} \right)}^{2} + {\delta\left( \tau_{r_{b}} \right)}^{2}}}} & (72) \end{matrix}$

In the next Section we will provide numbers for the relative error bounds we derived so far.

D. Measurements with Photons and Phonons

Let us see how well we can measure proper time, red shift and the Schwarzschild radius of the Earth using the electromagnetic field in a cavity or the phononic field in BEC. The change between the two systems can be readily obtained by substituting the speed of light c with the speed of sound c_(s0). The physical implementation, however, carries different challenges. This is due to the fact that both the amount of squeezing and the lifetime of the excitations are significantly different for photons and phonons within state-of-the-art experimental set-ups.

From Eq. (65), we see that the precision increases with the frequency ω_(0,k)=νπk/L_(p) of the mode under consideration. Hence, we choose a small cavity length L_(p)=200 μm. For photons, we have ν=c and we can assume a maximal squeezing of r=1.7 which is about 15 dB—the highest squeezing achieved to date. Using the reflectivity of current state-of-the-art mirrors, the largest time scale for photons is of the order of τ_(t,photon)=70 ns, which corresponds to the time that it takes to decrease the initial photon number n_(i) to n_(i/e). We can assume a power of 1 mW in the cavity, which corresponds to about 10⁴ photons for a wavelength of 500 nm. For the phonon field, we have v=c_(s0), which can be of the order 0.05 m/s, and we assume the maximal squeezing to be r=5.3 and the maximal number of phonons 10⁴. In a BEC, the maximum time scale is the decoherence time of the phonons of about τ_(t,photon)=6 s. With these numbers, we find for the photon field an error bound δ_(h,k) ^(SMSC)(τ_(r) _(b) ) of the order of 10⁻⁷/(√{square root over (M)}k) and for phonons δ_(h,k) ^(SMSC)(ρ_(r) _(b) )˜10⁻¹⁰/(√{square root over (M)}k). Since the same time interval cannot be measured twice, we can set M=1. We see that the phononic field performs better than the light field although the frequencies are much higher for the light field. One reason for this to occur is that it is not the frequency but the maximal achievable phase change ψ_(k)=ω_(0,k)τ that governs the relative error bound. The phonon states live longer than the photon states. But this would only lead to equivalence of the error bounds in our case. Another reason is that one can achieve much higher squeezing in a phononic field.

$\begin{matrix} {\mspace{79mu}{{{{\Delta_{h,k}^{SMSC}\left( {\tau\text{?}} \right)} = {\tau\text{?}\text{?}_{h,k}^{SMSC}\left( \tau_{r_{b}} \right)}}\mspace{20mu} 10^{- 16}s\text{/}\left( {\sqrt{M}k} \right)}\mspace{20mu}{10^{- 9}{s/\left( {\sqrt{M}k} \right)}}\mspace{20mu}{r_{b} ⪢ r_{S}}\mspace{20mu}{{2\frac{r_{b}}{r_{S}}{f\left( r_{b} \right)}{{1 - {f\left( r_{b} \right)}^{1/2}}}} \approx 1.}{\text{?}\text{indicates text missing or illegible when filed}}}} & (73) \end{matrix}$

We can conclude that the values of the relative error bounds for the estimation of the Schwarzschild radius and those for the estimation of the red shift are of the same order of magnitude. Therefore, we will only discuss the estimation of the Schwarzschild radius in the following.

Assuming r_(b)=6371 km and r_(s)=8.89 mm, we find the value 1.43×10⁹ for the factor (73) that would appear in Eq. (72), which has to be multiplied with the error bounds for the proper time measurements in order to obtain the relative error bounds δ(z_(∞)) for the estimation of the Schwarzschild radius. In the gravitational field of the Earth we have r_(b)>>r_(S) and we can assume that both clocks reach their maximal accuracy by approximately the same proper time at r_(∞). Considering the SMSC state and the mode k, we find a relative error bound δ_(k) ^(SMSC,∞,light)(r_(S)) of the order of 10³/(√{square root over (M)}k) for the light field and ϵ_(k) ^(SMSC,∞,BEC)(r_(s))˜1/(√{square root over (M)}k) for the phonon field. In contrast to the measurement of proper time, the estimation of the Schwarzschild radius can be performed several times. Let us assume that the experiment is performed for the lifetime τ_(l) of the states, which we assumed to be 70 ns for the light field and 6 s for the phononic field. If we assume that the experiment will be restarted without delay, we can set M=τ/τ_(l). We find for a measurement time of one year, i.e. τ=3.17×10⁷ s a relative error bound δ^(SMSC,∞,light)(r_(s)) of the order of 10⁻⁶/k for the light BE field_(s) in a cavity and a relative error bound δ_(k) ^(SMSC,∞,BEC)(r_(S)) of the order of 10⁻⁴/k for the phonon field in a BEC. This shows that the photonic clock is more precise for the measurement of the Schwarzschild radius than the phononic clock. In general, we find for the quotient of the two error bounds

$\begin{matrix} {\mspace{79mu}{{{\frac{\delta_{k}^{{SMSC},\text{?},{light}}\left( r_{S} \right)}{\delta_{k}^{{SMSC},\infty,{BEC}}\left( r_{S} \right)} \sim 10^{- 2}}\mspace{20mu} r\text{?}}\mspace{14mu}\mspace{20mu} r_{b}\mspace{20mu}\text{?}\mspace{14mu}\mspace{20mu} r_{b}{\text{?}\text{indicates text missing or illegible when filed}}}} & (74) \end{matrix}$

However, if we use two differently oriented extended clocks, we are able to let them interfere locally circumventing the transmission of any signal. We will discuss this set up in the next sections.

Section VI. Local Estimation of the Schwarzschild Radius

In this section, we investigate how to estimate the Schwarzschild radius locally by using two differently oriented clocks. First, we consider the comparison of independent readings of the clocks. At the end, we will also consider measurements of the whole system.

$\begin{matrix} {\mspace{79mu}{{{\psi\text{?}} = {\omega\text{?}\left( {1 + \Theta} \right)\tau\text{?}}}\mspace{14mu}\mspace{20mu}{{\psi\text{?}} = {\omega_{0,k}\tau_{r_{b}}}}\mspace{20mu} r_{b\mspace{11mu}}\mspace{20mu}{\tau_{eff} = {\left( {1 + \Theta} \right)\tau_{r_{b}}}}\mspace{20mu}\tau_{r_{b}}\mspace{20mu}\tau_{b}{{\delta^{local}\left( r_{S} \right)} = {{\frac{1}{r_{S}}{{\frac{\partial}{\partial r_{S}}\ln\frac{\tau_{eff}}{\tau_{r_{b}}}}}^{- 1}\sqrt{{\delta\left( \tau_{eff} \right)}^{2} + {\delta\left( \tau_{r_{b}} \right)}^{2}}} = {\frac{1}{r_{S}}{\frac{\partial\Theta}{\partial r_{S}}}^{- 1}\sqrt{{\delta\left( \tau_{eff} \right)}^{2} + {\delta\left( \tau_{r_{b}} \right)}^{2}}}}}{\text{?}\text{indicates text missing or illegible when filed}}}} & (75) \end{matrix}$

Considering r_(b)>>r_(s) and the SMSC state, it follows that

$\begin{matrix} {{\delta_{k}^{{cav},{loc}}\left( r_{S} \right)} \approx {\frac{4r_{b}^{2}}{r_{S}L_{p}}\frac{c_{rod}^{2}}{c^{2}}\sqrt{2}{\delta_{k}^{SMSC}\left( \tau_{r_{b}} \right)}}} & (76) \end{matrix}$

for the light field in the cavity and

$\begin{matrix} {{\delta_{k}^{{BEC},{loc}}\left( r_{S} \right)} \approx {\frac{4r_{b}^{2}}{r_{S}L_{p}}{{{\frac{3}{2}\frac{c^{2}}{c_{rod}^{2}}} - {\frac{c^{2}}{c_{s\; 0}^{2}}\frac{\delta_{p}^{\prime}L_{p}}{8}}}}^{- 1}\sqrt{2}{\delta_{k}^{SMSC}\left( \tau_{r_{b}} \right)}}} & (77) \end{matrix}$

for the phonon field in the BEC. δ′_(p)L_(P) is the relative change of the density of the BEC over the length of the trapping potential due to the difference between the external linear potential ΔV and the linear part of the gravitational field. We find that the sensitivity increases with increasing δ′_(p)L_(P). This is because δ′₀L_(p) appears quadratic in the relative frequency shift Θ. Let us assume that δ′_(p)L_(P)˜10⁻². Then, we obtain that c²δ′_(p)L_(P)/8c_(s0) ² is of the order 10¹⁸. As stated above, for an aluminium rod, c²/c² _(rod) is of the order 10⁹. Therefore, the effect of the density perturbation in the BEC dominates and we can write

$\begin{matrix} {{\delta_{k}^{{BEC},{loc}}\left( r_{S} \right)} \approx {\frac{4r_{b}^{2}}{r_{S}L_{p}}\frac{c_{s\; 0}^{2}}{c^{2}}\frac{8}{\delta_{p}^{\prime}L_{p}}\sqrt{2}{{\delta_{k}^{SMSC}\left( \tau_{r_{b}} \right)} \cdot \Delta}\; V}} & (78) \end{matrix}$

Employing two extended clocks at the same position can have an additional advantage. They can be brought into interference. This means that one could entangle the modes before the measurement. Furthermore, quantum mechanical measurements could be performed on both systems without the necessity of encoding the state of one of the two clocks into an intermediate quantum system, such as a laser beam, for the transportation of the quantum state. We will discuss this possibility in the next section by giving an explicit measurement that can approach the error bound given by the quantum Fisher information

Section VII. A Phononic Mach-Zehnder Interferometer

Mach-Zehnder interferometers on BECs have been experimentally proposed in early work such as T. Berrada, S. van Frank, R. Bucker, T. Schumm, J. F. Schaff, and J. Schmiedmayer, Nature Communications 4, 2077 EP (2013). and J. Grond, U. Hohenester, J. Schmiedmayer, and A. Smerzi, Phys. Rev. A 84, 023619 (2011, when the condensate is trapped in a harmonic potential. Although experimental interferometer schemes with box-trapped phonons as the quantum states are yet to be developed, proposals such as A. C. J. Wade, J. F. Sherson, and K. Mølmer, Phys. Rev. Lett. 115, 060401 (2015) that allow the mapping of the quantum state of a BEC into a light field make this idea feasible in the near future.

-   -   ω_(h,k)

ω_(r,k)=ω_(h,k)(1+Θ)

ω_(v,k)−ω_(b,k)=ω_(h,k)Θ

The full transformation of the state depicted in FIG. 4 is given by

$\begin{matrix} {\left. \;{{\hat{B}\text{?}\left( {{\hat{U}\left( {\phi^{\Delta}\left\lbrack {\Theta\left( r_{S} \right)} \right\rbrack} \right)} \otimes {\hat{U}\left( {\pi\text{/}2} \right)}} \right)}❘\Psi_{0}} \right\rangle,} & (79) \\ {\left. {{{\left. {\left. {\left. {{{{\left. {\left. \mspace{79mu}{{\pi\text{/}2}❘\Psi_{0}} \right\rangle = {B\text{?}{\Psi_{h}}\Psi\text{?}}} \right\rangle\text{?}^{\Delta}{{\Theta\left( r_{S} \right)}}} = {{{\psi\text{?}} - \psi_{k}} = {\Theta\text{?}}}}\mspace{20mu}{{B\text{?}} = {{\text{?}\mspace{14mu}\text{?}\text{?}} = {\frac{1}{2}\left( {{\hat{a}\text{?}{\hat{a}}_{h}} - {a\text{?}\text{?}}} \right)}}}\mspace{20mu}\left. a_{h}\rightarrow{\text{?}\text{?}\text{?}\text{?}} \right.}\mspace{14mu}\mspace{20mu}{\left. {\hat{a}}_{r}\rightarrow{\hat{a}}_{v} \right.\mspace{20mu} ❘\Psi_{h}}} \right\rangle = {{{\hat{D}(a)}\left. 0 \right\rangle} = {❘a}}} \right\rangle\mspace{11mu}\mspace{20mu} ❘\Psi_{v}} \right\rangle = {{{S(r)}\left. 0 \right\rangle} =}}}r} \right\rangle\mspace{20mu}{{II}:={\hat{B}\text{?}\overset{.}{M}}}\mspace{79mu}{\overset{.}{M} = {{\hat{J}\text{?}} = {\frac{1}{2}\left( {{\hat{a}\text{?}{\hat{a}}_{h}} - {\hat{a}\text{?}{\overset{.}{a}}_{v}}} \right)}}}\mspace{20mu}{{{\delta_{\exp}^{MZI}\left( r_{S} \right)} \approx {\frac{1}{\psi_{0,k}\sqrt{M}}\frac{4r_{b}^{2}}{r_{S}L_{p}}\frac{c_{\Delta\bigcap}^{2}}{c^{2}}\frac{8}{\delta_{p}^{\prime}L_{p}}{\Delta_{\phi}}^{\Delta}}},\mspace{20mu}{\Delta_{\phi}}^{\Delta}}} & (80) \\ {\mspace{79mu}{{{\Delta\phi}^{\Delta} = \sqrt{\frac{{{\alpha }^{2}c^{{- 2}r}} + {\sinh^{2}(r)} + {{\tan^{2}\left( \phi^{\Delta} \right)}{F\left( {\alpha,r} \right)}}}{\left( {{\alpha }^{2} - {\sinh^{2}(r)}} \right)^{2}}}}\mspace{20mu}{{F\left( {\alpha,r} \right)} = {{\alpha }^{2} + {2{\sinh^{2}(r)}\text{?}^{2}(r)}}}\mspace{20mu}{\alpha = {i{\alpha }}}\mspace{20mu}{\text{?}/2}\mspace{20mu}{\phi_{opt}^{\Delta}\hat{=}0}\mspace{20mu}{{\alpha }^{2}\text{?}{\sinh^{2}(r)}}}} & (81) \\ {\mspace{79mu}{{{\delta_{{QF}\; 1}^{MZI}\left( r_{S} \right)} \approx {\frac{1}{\psi_{0,k}\sqrt{M}}\frac{4r_{b}^{2}}{r_{S}L_{p}}\frac{c_{s\; 0}^{2}}{c^{2}}\frac{\delta}{\delta_{p}^{\prime}L_{p}}{\Delta\phi}_{{QF}\; 1}^{\Delta}}},{\text{?}\text{indicates text missing or illegible when filed}}}} & (82) \end{matrix}$

where (as per L. Pezze et all, above)

$\begin{matrix} {{{{\Delta\phi}_{{QF}\; 1}^{\Delta} = {{\frac{1}{\sqrt{{\alpha }^{2} \in^{2}{{+ \sin}\mspace{11mu}{h^{2}(r)}}}} \cdot {\alpha }^{2}} \sim {\sinh^{2}(r)}}}2{{{esc}\left( {2\theta} \right)}}\sqrt{N_{ph}/N_{nt}}}\mspace{14mu}\theta} & (83) \end{matrix}$

A. Comparison with the Optimal Bound Given by the QFI

It is not possible to use the exact same parameters to compare the sensitivities achieved with phonon number difference and with parity/QFI since each reach their minimum with different values of the squeezing parameter. For the Ĵ_(z), measurement, a value of |α|=10³ (corresponding to a number of Rubidium atoms of 10⁹) fixes the optimal sensitivity with a squeezing parameter of r=3.77. Using equation (80) we obtain a relative error of δ_(exp) ^(MZI)(r_(S))=4.92×10⁻⁸ for a single-shot and δ_(exp) ^(MZI)(r_(S))=2.14×10⁻¹¹ after one year of integration time. For the parity measurement/QFI the relative error reads δ_(QFI) ^(MZI)(r_(S))=3.43×10⁻⁸ after a single shot-measurement, further leading to δ_(QFI) ^(MZI)(r_(S))=1.49×10⁻¹¹ after a year of integration, with a squeezing parameter of r=5.29. The additional correction assuming a pumped-up SU(1,1) setup was chosen with a mixing angle of θ˜0.1, which pumps the number of phonons up to 1% the number of atoms.

B. Comparison with Previous Schemes

One of the main advantage of our proposal is its dimensions. At the scale of micrometres, its portability exceeds that of current atom interferometry setups by at least two orders of magnitude.

While its sensitivity appears be below than those reported by the other proposals, it is also important to note that the time scale at which the experiment can be realised is shorter than their integration time, since the reported running time of our proposal corresponds to a single run. We can calculate the integration time T_(int)=M×6 s in order to reach the other setups' sensitivities, obtaining T_(int)≈1.96 hours to reach the precision of the atom interferometers and T_(int)≈3.22 hours to reach the precision of a previous on-chip BEC. As an additional note, another promising setup of length 10⁻⁵ m that relies on optomechanics was recently proposed. Although its relative error is yet to be experimentally confirmed, it is reported to be 10⁻¹³ in a cycle time of 10⁻³ s. However, our proposal is a purely static BEC, which further simplifies its experimental implementation.

Section VIII. Conclusions and Discussion

We have discussed above how to employ the phonon field in an effectively 1-dimensional BEC as a quantum probe of the gravitational field of the Earth. We compared it to the light field in a 1-dimensional optical cavity used for the same purpose. The external field can be easily tuned by finding the point at which the BEC bulk keeps its position.

Each single mode of the quantum fields can be seen as a spatially extended clock. We found that the BEC clock and the optical cavity clock can measure time with a relative precision of up to 10⁻¹⁶ and 10⁻⁷ for a single measurement, respectively. However, to achieve these levels of sensitivity, the duration of a single measurement on the BEC clock must be 6 s, the maximal lifetime of the phonons in the BEC, while a single measurement on the light field in the cavity only needs to take 70 ns, the life-time of the photon state in the cavity. In particular, the maximal time resolution of a light cavity clock is 10⁻¹⁶ s while it is about 10⁻⁶ s for the BEC clock. The reason for this is the much higher frequency of the light field modes in comparison to the phonon field modes. We assumed optical frequencies for the light cavity. Therefore, it is a proper lower bound for realistic detector systems. We also discussed the estimation of red shift and the Schwarzschild radius using two independent time measurements at different radial distances from the surface of the Earth. For an integration time of one year, we found a relative error bound of 10⁻⁶ and 10⁻⁴ for the light cavity clock and the BEC clock respectively.

We introduced a novel idea of how to measure the proper acceleration or the Schwarzschild radius of the Earth with two differently oriented BECs/cavities on the length scale of 200 μm. The basis of the measurement is the difference of the fundamental frequency of phonon modes in differently oriented BECs (or photon modes in differently oriented cavities) due to the deformation of the length of the trapping potential of the BEC (or the length of the cavity). There are three separate effects that play a role. The first is the fundamental difference between the proper length of a rigid rod and the radar length, which defines the fundamental frequency in a curved spacetime. This effect is completely given by the finite extension of the clocks we consider. The second effect is the deformation of the material which holds the cavity and the material that holds the devices that generate the external potential which holds the BEC, respectively. Finally, the third and dominant effect is a possible deformation of the BEC cloud due to a mismatch on gravitational parameters. We assumed the material to be a rod of mass density ρ_(rod) and Young's module E, and we assumed the whole device to be supported from below. We found that the effect of the deformation of matter is much larger than the effect of the deformation of spacetime for any realistic situation.

The measurement for the estimation of the gravitational acceleration or the Schwarzschild radius can be done by two independent time measurements on the two clocks and subsequent comparison of the results or by a measurement on the combined system of the two differently oriented clocks in a Mach-Zehnder interferometer-like scheme. The relative error bound for the estimation of the Schwarzschild radius on the length scale of 200 μm turned out to be, respectively, 10⁻⁵ and 10⁻¹⁶ for the light cavity in the optical regime and for the BEC if the rod which keeps the length of the BEC/cavity fixed is made from aluminium. Finally, for the MZI setup we get bounds of 10⁻¹¹ both for the Ĵ_(z) measurement and for the QFI/parity measurement with the additional improvement of pumped-up setups.

Detailed Discussion: Gradiometer

The collective oscillations of the atoms in BECs can also/alternatively be used to estimate a field gradient causing perturbations in a BEC, for example a gravity gradient. In principle, this can enable the measurement of a gravity gradient due to a mass. Below we demonstrate our method by considering spherical masses ranging from the size of the earth to small spheres with masses down to hundreds of milligrams. It will be appreciated that the method can be adapted for use with other masses of different sizes and shapes. It is notable however that the measurement scheme functions on the length scale of a BEC, which is of the order of micrometres.

The principle of gradiometry is based on the fact that an external force adds a quadratic term to the trap potential which changes the frequencies of the collective oscillations of the BEC. Consider a spherically symmetric trapping potential given by V_(trap)=mω₀ ²(x²+y²+z²)/2), where ω₀ is the is the trapping frequency and m is the mass of the atoms in the BEC. In the case of the gravitational field of a massive sphere, the quadratic term in the potential is proportional to the gravity gradient. The BEC 38 is therefore spherically symmetric, and possesses a radius

₀.

If this setup is placed in a gravitational field represented by a Newtonian potential ϕ, the total potential is given by V=V_(trap)+mΦ. The gravitational field of a spherically symmetrical, homogeneous mass distribution can be written as Φ=−MG/R where M is the total mass, G is Newton's gravitational constant and R is the radial distance from the centre of the mass distribution. Assuming that the centre of the mass distribution is located along the z-axis at a distance R from the centre of the trap potential, the gravitational field can be written as Φ=−MG(z/R²−(ρ²−2z²)2R³) up to second order in the spatial distances, where ρ³=x²+y². Therefore, the total potential seen by the BEC can be written as V=m(ω_(⊥) ²ρ²+ω_(z) ²(z−z_(g))²)/2+C, where ω_(⊥)=(ω₀ ²+MG/R³)^(1/2), ω_(s)=(ω₀ ²−2MG/R³)^(1/2), z_(g)=MG/R²ω_(z) ² and C=−m(MG/R²)²/2ω_(z) ². The linear part of the Newtonian potential leads to a shift of the equilibrium position of the BEC.

By defining z′=z−z_(g) and by neglecting the constant term C (or cancelling it by an additional global contribution to the time evolution of phonon modes), a new potential is obtained:

$\begin{matrix} {V = {\frac{m}{2}{\left( {{\omega_{\bot}^{2}\rho^{2}} + {\omega_{z}^{2}z^{\prime 2}}} \right).}}} & (1) \end{matrix}$

It can be seen that the trap frequencies depend on the gravity gradient ∈_(grad)=2MG/R³, which forms the basis of its measurement.

In an axially symmetric potential such as Eq. 1, the BEC's stationary density can be approximated as n=μ(1−(ρ²+λ²z′²)/

²)U₀, where λ=ω_(a)/ω_(⊥) encodes the shape of the BEC (in our case λ<1 which corresponds to a prolate shape),

is the radius of the BEC, U₀=4πR²a_(scatt)/m, a_(scatt), is the scattering length and μ is the chemical potential. The radius of the BEC can be given in terms of the chemical potential as

²=2μ/mω_(⊥) ² and the total number of atoms can be derived from the density as N_(o)=8π

³μ/15U₀λ. This leads to an expression for the radius of the BEC as

=(15U₀λ)/4πmN_(a)ω_(⊥) ²)^(1/3).

The above approximation is called the Thomas-Fermi approximation and it describes the density profiles in experimental situations quite accurately if N_(a)a_(scatt)/ā_(O)>>1, where N_(a) is the number of atoms in the BEC, ā_(O)=(h/mω)^(1/2) and ω=(ω_(⊥) ²ω_(s))^(1/3). The parabolic density profile and the spheroidal shape of the BEC give rise to a specific set of modes of density perturbations the density perturbations have the spatial dependence δn∝r^(l)Y_(lm)(θ,ϕ) in spherical coordinates r=ρ²+t², θ=arccos(z/r) and ϕ=arctan(y/x)), where Y_(lm)(θ, ϕ) is the spherical harmonic function with angular momentum l and the z-component of the angular momentum vector m=±l or m=±(l−1). These two types of modes are denoted as δn_(l,l) and δn_(l,l−1). The corresponding angular frequencies are given as ω_(l,l) ²=lω_(⊥) ² and ω_(l,l−1) ²=(l−1)ω_(⊥) ²+ω_(s) ², respectively. The gravity gradient is expected to be much smaller than that of the trap potential. Therefore, the frequencies can be approximated as ω_(l,l)≈√{square root over (l)}(ω₀+GM/2R³ω₀) and ω_(l,l−1)≈l(l−3)GM/2lR³ω₀).

δn_(l,l) δn_(l,l−1) contain information about the gravity gradient ∈g_(rad)=2MG/R² as an external parameter. Therefore, it is possible to estimate its value from a measurement on these modes.

A lower bound for the relative error δ_(e) of any estimation of a given parameter ∈ imprinted on a given mode can be obtained from the quantum Cramer-Rao bound (QCRB). The QCRB can be expressed in terms of the quantum Fisher information (QFI) H_(e) and the number of measurements N as Δ_(e)=1/√{square root over (NH_(e))}. The corresponding relative error bound is obtained as δ_(e)=1/|∈|√{square root over (NH_(e))}. Since MG/R³ is imprinted only on the frequency, we obtain the minimal relative error for the estimation of ∈_(grad) by Gaussian error propagation from the minimal error of an estimation of a phase change Δϕ=Δωl as δ_(r)=δ_(Δϕ)/|d ln(Δϕ)/d ln ∈_(grad)|, where Δω=√{square root over (l)}∈_(grad)/4ω₀ for the mode δn_(l,l−1). Therefore, we find that δ_(c)=δ_(Δϕ). If the initial state of the mode under consideration is assumed to be a Gaussian state and if only a one single mode is being measured, the optimal precision for the measurement of a phase change is reached for a squeezed vacuum state and the corresponding QFI becomes H_(Δϕ)=δmin+1), where n=sinh²r is the number of squeezed phonons, where r is the squeezing parameter. Then:

$\begin{matrix} {\mspace{79mu}{{{{\delta\text{?}} = {{\frac{2\omega_{0}^{2}}{{\alpha(l)} \in_{grad}}{\frac{1}{\sqrt{l}\omega_{0}t\sqrt{2{{Nn}\left( {n + 1} \right)}}} \cdot {\alpha(l)}}} = 1}}\mspace{85mu}{\delta\; n_{l,l}}\mspace{79mu}{{\alpha(l)} = {{\left( {l - 3} \right)\text{/}l}}}}\mspace{79mu}{\delta\; n_{l,{l - 1}}}{\text{?}\text{indicates text missing or illegible when filed}}}} & (2) \end{matrix}$

We now evaluate the relative error bound in Eq. (2) for two specific situations; the gravitational field of the earth and the gravitational field of a sphere of tungsten or gold. Lets assume that the trapping frequency is a ω₀=2π×0.2 Hz. For a rubidium-87 BEC of 10⁶ atoms this would lead to a radius of 120 μm and a central density of n n(0)˜10¹¹ cm⁻³, which is fully in the currently experimentally accessible regime. This would lead to a radius of 380 μm and a central density of n(0)˜10¹² cm⁻¹. The small density of the BEC is an advantage as it leads to a long half life time of the BEC density. It can be was shown that the density depends on time as dρ(t)/dt=˜Dρ(t)², where D is the decay constant. Therefore, after solving the differential equation, we find a quadratic dependence of the density half life time on the inverse density, i.e. t_(l,l)=3/2Dρ². Previous experiments with rubidium atoms found the decay constant to be D=1.8×10⁻²⁰ cm⁰s⁻¹. For a density of the order of 10¹² cm⁻³, this leads to a theoretical half life time of the order of the order of 10⁵ s. Therefore, it can safely assumed that the BEC density can be kept constant for a duration of each experiment of the order of 100 s. Then, the second limiting time scale that one has to consider is the coherence time of the phonons. For BECs of temperatures T below or of the order of the chemical potential divided by the Boltzmann constant k_(B), the damping rate is given by γ˜√{square root over (l)}ω₀(k₀T/μ)^(3/2)(n(0)a_(scatt) ³)^(1/2). Assuming l=3 and a temperature of the BEC of 0.5 nK, which can be achieved in experiments, the inverse damping rates are of the order of 10² s and 10³ s for 10⁶ and 10⁸ rubidium atoms, respectively. From the inverse damping rate, one can calculate the coherence time of a squeezed phonon state. For the above parameters, the coherence times are calculated to be between 10² m and 10³ m seconds. Therefore, it is justified to assume a duration of 100 s for each experiment.

For the number of independent consecutive measurements, a value of 10⁴ is assumed which corresponds to a total measurement time of about one and a half weeks. The number of squeezed phonons in the cases of 10⁶ and 10⁸ atoms, respectively, is set to n=10³ and n=10⁴ (which corresponds to a squeezing parameter of r˜4 and r˜5, respectively). The mode number l=3 is used. For the case of the gravitational field of the earth and R of the same order as the radius of the earth, it is found that ω₀ ²/∈_(grad) is of the order of 10⁶. This leads to a relative error bound of the order of 10⁻² and 10⁻³ for 10⁶ and 10⁸ atoms, respectively. Hence in principle, it is possible to measure the gravitational gradient due to the gravitational field of the earth on the length scale of a BEC using the phonons in the condensate.

A 100 mg gold or tungsten sphere has a radius of the order of 1 mm, and one can assume a distance between the centre of the BEC and the centre of the sphere of the same order. With ω₀=2π×0.2 Hz, it is found that ω₀ ²/∈_(grad) is again of the order 10⁶ and a relative error bound of the order of 10⁻² and 10⁻³ for 10⁶ and 10⁸ atoms, respectively, by using the same system parameters as above. The reason for this scaling is that by considering R to be always of the same order as the radius of the sphere independently of its mass, the expression 2MG/R³˜2ρG is obtained, where ρ is the mass density of the sphere.

Therefore, the absolute sensitivity of the scheme is given as ∈_(grad)10⁻²˜10⁻⁷ s⁻² and ∈_(grad)10⁻³˜10⁻⁸ s⁻² for 10⁶ and 10⁸ atoms in the BEC, respectively. On the length scale of 100 μm, such small gradients correspond to gravitational field differences of 10⁻⁹ Gal (10⁻¹² g) and 10⁻¹⁰ Gal (10⁻¹³ g), respectively. 10⁴ repetitions of the experiments are considered. Therefore, the single shot sensitivity would be comparable to a differential force sensitivity of the order of 10⁻⁷ Gal (10⁻¹⁰ g) and 10⁻⁸ (10⁻¹¹ g), respectively, which is one order of magnitude better than what is achieved by state of the art atom interferometry.

To measure the gravity gradient induced by even smaller masses, the trapping frequency would have to be increased to decrease the radius of the BEC. This would lead to a decrease of the sensitivity. For example, a gold or tungsten sphere with a mass of the order of 100 μg has a radius of the order of 100 μm. To get a BEC with a radius of the order of 10 μm with 10⁵ rubidium-87 atoms, we have to use a trapping frequency ω₀=2π×20 Hz. The corresponding central density is of the order of 10⁻¹³ cm, which still allows for experimental times of 100 s in principle. Using the same parameters for the phonon state as above, we obtain a relative error bound of the order of 10⁻¹ and 10² for 10⁶ and 10⁸ atoms, respectively and an absolute error bound of 10.8 gal (10⁻¹¹ g) for 10⁶ atoms and 10⁻⁹ gal (10⁻¹² g) for 10⁸ atoms with 10⁴ repetitions.

To have a chance to reach the relative error bound given in Eq. (2), there has to be a phase reference to compare the phase change of the probe state with. In particular, the frequency of the reference has to be exactly √{square root over (l)}ω₀. One possibility to avoid this would be to use two modes of perturbations of the BEC density and compare them with each other such as in a Mach-Zehnder interferometer (MZI) scheme. By comparing two modes with the same angular momentum l but different quantum number m, the difference of the two frequencies can be approximated as:

$\begin{matrix} {{\Delta\omega}_{i} = {{\omega_{l,l} - \omega_{l,{l - 1}}} ≉ {\frac{3 \in_{grad}}{2\sqrt{l}\omega_{0}}.}}} & (3) \end{matrix}$

The frequency difference would lead to an accumulated phase difference between the modes Δϕ_(l)=Δω_(l)t. Note that the frequency difference (Eq. 3) decreases with increasing angular momentum l. This is in contrast to the change of frequency of each single mode due to the gravitational field considered above. This would enable multi-mode operations that can be used for the implementation of MZI schemes. In metrology with optical modes, there are two distinct classes of schemes of Mach-Zehnder interferometry commonly used; SU(2) and SU(1,1). SU(2) schemes only consider passive optical elements and the SU(1,1) schemes contain active optical elements. In all schemes the optimal QFI depends on the total number of particles in the two arms of the interferometer which we denote as n in the following. It can be shown that the optimal QFI for a phase measurement for the SU(2) scheme is reached for two identically squeezed and coherently displaced states at the two input ports of the first beam splitter of the MZI. The total number of squeezed particles has to be ⅔ of the total number particles n. For the case of n>>1, the optimum becomes H_(δϕ) ^(S(1)(2))≈8n(n+2)/3. In contrast to the SU(2) scheme, the SU(1,1) scheme uses optical parametric amplifiers (OPA) instead of passive beam splitters. At the OPA the light fields interact non-linearly and they are additionally squeezed. It can also be shown that the optimal QFI for the SU(1,1) scheme is given as H_(δϕ) ^(SU(1,1))≈4{circumflex over (n)}({circumflex over (n)}+2)/3 for n>>1. To reach this QFI, two coherent states with the same number of particles have to injected into the two ports of the MZI and the number of squeezed particles that are created at the OPA should be ⅔ of the number total number of particles in the arms of the interferometer. If the number of squeezed particles is fixed, H_(δϕ) ^(SU(2)) and H_(δϕ) ^(SU(1,1)) differ from the optimal single mode QFI that we discussed above approximately by a factor ¾ and ⅜, respectively.

However, it is in general more easy to create highly excited coherent states than large squeezing. A particular adaptation of the SU(1,1) scheme is also beneficial when the number of squeezed particles is much smaller than the total number of particles. It is called “pumped-up” SU(1,1) interferometry because the active optical elements are additionally pumped by a strong coherent light field. It is shown that the optimal QFI for the pumped up SU(1,1) scheme becomes H_(δϕ) ^(pu)=ne^(2r)/4. For n=sinh²r and n>>1, and it is found that differs from the one mode QFI by a factor 1/16. However, in situations where the number of coherent particles is much larger than the number of squeezed particles, the application of the pumped up SU(1,1) interferometry can be truly beneficial.

SUMMARY

The above work shows that the density of a trapped BEC is modified by gravity and that this affects the evolution of collective BEC modes. In the case of trapping potentials that produce a constant and uniform density in the absence of gravity (e.g. box potentials), one can find modes that acquire different phase evolution in the presence of gravity.

A phase acquired by a quantum state |φ

during its evolution is given by e^(iθ(t))|φ

. In this case the phase shift θ(t) is a function of time t. In our work, the quantum state can describe one or more modes. In the discussions above we focus on the state of two modes but our scheme can be applied to more modes if required.

To measure a phase we need it to be a relative phase i.e. of the form

|φ_(f)

=|φ₁

+e ^(iθ(t))|φ₂

Where |φ₁

and |φ₂

are two different states of the system.

As explained above, the phase shift θ(t) depends on the gravitational parameters and that is why we can then us it to measure gravity and its gradient.

It will be appreciated that the case described above is a special one. As well as affecting the phase difference between two modes of a system, an external field applied to a closed system might cause excitations to transfer between two modes of the system, or new excitations to be created within one or more of the modes.

A more general evolution of a quantum system can be modelled using a unitary transformation. A unitary is given by e^(i{circumflex over (Θ)}(t)) where now {circumflex over (Θ)}(t) is not only a function of t but a matrix. This is the most general evolution a quantum system undergoes when the system is closed (when it does not interact with the environment). It acts on the state of the system as e^(i{circumflex over (Θ)}(t))|φ

where the quantum state |φ

can describe one or more modes.

|φ_(f)

=e ^(i{circumflex over (Θ)}(t)) |φ

=ΣC(t)_(ij)|φ_(i)

|φ_(j)

Here C(t)_(ij) are functions of t. A phase is then a special case where C(t)_(ij)=0 if i≠j and C(t)_(ij)=θ_(i)(t) if i=j. Here θ_(i)(t) is again a function of t. If the modes are phonon modes trapped in a potential in the presence of gravity, the numbers C(t)_(ij) will depend on the gravitational parameters and can be used to measure them. In M. Ahmadi, D. E. Bruschi, C. Sabin, G. Adesso and I. Fuentes, Sci. Rep. 4, 4996 (2014) (noted in the background section above) we analyse the case of an oscillating cavity in flat spacetime with acceleration a.

If the unitary transformation acts on a m number of modes but we measure only n (with n<m) or if the system undergoes losses due to the interaction with the environment, then the modes of interest n evolve under a Complete Positive map which is the transformation we obtain after tracing out the degrees of freedom which we do not measure.

In both cases, the evolution of the modes is no longer pure (pure means it can be represented by a vector |φ

in the Hilbert space). Instead the states are mixed, which means that the state can no longer be represented by a vector in the Hilbert space but by a matrix, called the density matrix, which is an operator acting on the Hilbert space. The schemes described herein work with both pure and mixed states.

For more general trapping potentials that produce time varying and/or non-uniform density, the mode evolution may be a complete positive map that encodes the gravitational parameters. Complete positive maps include the case where the evolution is a phase (as discussed above), but also extend to unitary transformations such as a Bogoliubov transformations and to non-unitary evolutions due to decoherence. The field parameters (including the gravitational field strength (gravimeter example above) and its gradient (gradiometer example above)) can be estimated in the Newtonian or General relativistic regimes by measuring the system (for example the phase difference between the modes) using absorption imaging, interaction with a cavity mode, interaction with a laser, interaction with impurities among other methods.

The scheme therefore considers two or more frequency modes (including zero momentum and sharp momentum modes, also superpositions of sharp momentum modes) of a single or multiple trapped BECs (trapping potentials can give rise to uniform, non-uniform, constant and varying density). The modes can be prepared in a suitable quantum state (which can be separate or entangled states and can include coherent, squeezed, or any other pure or mixed state) such that the mode evolution encodes the field (e.g. gravitational) parameter to be estimated.

We thus disclose a method of measuring the gravitational field using one or more trapped BECs and at least two atomic field modes within those BECs, wherein the gravitational field leads to a change of the field modes' time evolution amounting to a gravitational-dependent phase or/and a number of excitations being exchanged between the modes (the number of excitations exchanged by the modes depending on the gravitational field) or/and a number of new excitations being created in the modes (the number of created excitations depending on the gravitational field). The method comprises:

-   -   selecting at least two modes whose initial states have changed         due to the gravitational field;     -   measuring the phase difference, or counting the number of         excitations exchanged by the modes or counting the number of         excitations created within the modes; and     -   using either the measured phase difference, the number of         excitations exchanged, the number of excitations created, or a         combination of any of these to infer the acceleration.

Number of excitations may be counted using standard methods including absorption imaging, homodyne or heterodyne detection, through interactions with a laser field or cavity field.

The above work differs from atom interferometry in that we consider the atoms to be trapped, and we include interactions between those trapped atoms so that the modes we consider are collective modes. Our modes are sharp in frequency and delocalised in the extent of the system. Instead atom interferometry considers free falling spatially separated single particle states.

The above work is distinguished from double-well systems in that such schemes consider only right atoms vs left atoms rather than phononic modes, and do not consider how gravity affects the density. 

1. A method of measuring an acceleration due to a gravitational field using one or more BECs, wherein the one or more BECs include: a first mode propagating in a first direction; and a second mode propagating in a second direction different to the first direction; wherein the method comprises: measuring a phase difference between the first mode and the second mode; and using the measured phase difference to measure the acceleration due to the gravitational field.
 2. The method of claim 1, wherein the first direction is parallel to a direction of the acceleration due to the gravitational field and the second direction is perpendicular to the direction of the acceleration due to the gravitational field.
 3. A method of measuring an acceleration using one or more trapped BECs and at least two atomic field modes within those BECs, wherein a density distribution of at least one of the one or more BECs is modified by the acceleration, which leads to a change of the field modes' time evolution; wherein the method comprises: selecting two modes that are differently affected by the acceleration; and measuring the acceleration-induced difference; and using the measured acceleration-induced difference to infer the acceleration.
 4. The method of claim 3, wherein the acceleration is an acceleration due to gravity, and wherein the acceleration-induced difference between the two modes comprises one or more of: (i) an acceleration-dependent phase difference, (ii) a number of excitations exchanged between the modes, and (iii) a number of new excitations created within one or both of the modes
 5. The method of claim 1, or claim 2 or claim 4, wherein the phase difference between the first mode and the second mode is used to estimate the local Schwarzschild radius r_(s).
 6. The method of claim 5, wherein the local gravitational field is estimated from the local Schwarzschild radius r_(s) using g=r _(s) c ²/2 r _(b) ² where r _(b) is the radial height coordinate within the gravitational field. condensing the at least one BEC in a trapping potential which is tuned to partially cancel the acceleration caused by the gravitational field, such that a density perturbation δ_(ρ) is created in the BEC.
 8. The method of claim 7, wherein the trapping potential is ${\Delta\; V} \approx {{- \frac{m^{2}c^{2}}{h^{2}}}\frac{r_{S}x^{\prime}}{r_{b}^{2}}}$ where  x^(′) = x− $\delta_{p} = {{- \frac{1}{2}}\frac{c^{2}}{c_{s\; 0}^{2}}\left( {\frac{r_{S}}{{\overset{\_}{r}}_{b}^{2}} - \frac{r_{S,0}}{{\overset{\_}{r}}_{b,0}^{2}}} \right)\left( {x^{\prime} - \frac{x_{t}^{\prime}}{2}} \right)}$ condensing a first BEC in a first atom trap having a first length L_(p) between the first and second locations (which may be boundaries) r _(b) and r _(t), condensing a second BEC in a second atom trap having a second length L_(p) between the third and fourth locations (which may be boundaries) z_(l) and z_(r), the second length being equal to the first length, wherein the first mode is measured in the first BEC and the second mode is measured in the second BEC.
 10. The method of claim 9, wherein the first BEC is confined in a box potential and wherein the second BEC is confined in a box potential.
 11. The method of claim 9 or claim 10, wherein the length L_(p) of each BEC is held constant by a respective rigid rod, and wherein deformation of the rigid rod due to gravity is taken into account for the first BEC but neglected for the second BEC.
 12. The method of any one of claims 8 to 11, further comprising creating a squeezed state for a mode k of the first BEC.
 13. The method of claim 12, further comprising preparing a coherent state or a squeezed state for a mode of the second BEC having the same mode number k as the squeezed mode of the first BEC.
 14. The method of claim 12 or claim 13, further comprising entangling the first and second modes.
 15. The method of claim 14, wherein a phase shift is introduced between the first and second modes, wherein the phase shift is optionally π/2.
 16. The method of any preceding method claim, wherein the first mode is a phonon mode.
 17. The method of any preceding method claim, wherein the spacetime of the or each BEC is described using Schwarzschild geometry in order to calculate the acceleration due to gravity.
 18. A quantum accelerometer operable to perform the method of any one of claims 1 to
 17. 19. A method of measuring a field gradient using one or more trapped BECs, each having a density distribution, and at least two atomic modes within those BECs; wherein the density distributions of the one or more BECs are modified by the field, which leads to a change of the atomic modes' time evolution; wherein the method comprises: selecting two atomic modes that are differently affected by the field; and measuring a field-induced difference between the selected atomic modes; and using the measured field-induced difference to infer a gradient of the field.
 20. The method of claim 19, wherein the field is a gravitational field, and wherein the field-induced difference between the two modes comprises one or more of: (i) a field-dependent phase difference, (ii) a number of excitations exchanged between the modes, and (iii) a number of new excitations created within one or both of the modes.
 21. The method of claim 19 or claim 20, wherein the method comprises using a single trapped BEC having a density distribution that is perturbed by the external field, and wherein the two differing modes are selected from modes existing in the single BEC.
 22. The method of any one of claims 19 to 21, wherein each atomic mode is defined by an angular momentum quantum number and a magnetic quantum number, and wherein measuring the field-induced difference between two modes involves selecting two modes with the same angular momentum quantum number, but different magnetic quantum numbers.
 23. The method of claim 22, wherein the magnetic quantum numbers of the selected modes differ by one.
 24. The method of any preceding method claim, wherein the or each BEC is trapped in a trapping potential.
 25. The method of claim 24, wherein the trapping potential comprises an electromagnetic portion and a field-induced portion, and wherein the field-induced portion is due to a gravitational field caused by a mass distribution in the vicinity of the BEC causing a distortion in the shape of the BEC.
 26. The method of claim 25, wherein the electromagnetic portion of the trapping potential is a harmonic trapping potential.
 27. The method of claim 26, wherein the gravitational field adds a linear term to the electromagnetic portion of the trapping potential.
 28. The method of any preceding method claim, wherein measuring the field-induced difference comprises taking a plurality of independent consecutive measurements of a phase difference between the modes, and wherein the measurements are averaged to produce an average phase difference.
 29. The method of any preceding method claim, where in one or more atomic modes within the BEC are prepared in a squeezed state.
 30. A quantum gradiometer operable to perform the method of any one of claims 19 to
 29. 31. A gradiometer for measuring a field gradient; the gradiometer comprising: a trapped BEC having a density distribution, wherein the trapped BEC comprises a plurality of atomic modes, and wherein the density distribution of the trapped BEC is modified by the field; a measurement system for measuring a field-induced difference between two selected atomic modes within the trapped BEC; and a processor for inferring a field gradient from the measured field-induced difference.
 32. The gradiometer of claim 13, wherein the field is a gravitational field, and wherein the field-induced difference between the two modes comprises one or more of: (i) a field-dependent phase difference, (ii) a number of excitations exchanged between the modes, and (iii) a number of new excitations created within one or both of the modes. 